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COMPUTER PROGRAM FOR THERMAL ANALYSIS OF 
TANK-MOUNTED MULTILAYER INSULATION 
by Glenn R. Cowgill and William R. Johnson 
Lewis Research Center 

SUMMARY 

A computer program was developed to analyze the space-hold thermal performance 
of various multilayer insulation (MLI) configurations on the upper half of an oblate 
spheroidal liquid-hydrogen tank. The thermal model considered was an unshrouded tank 
within a hypothetical sun-oriented vehicle such that the payload was the only heat source 
in the system. 

The thermal code CINDA-3G (Chrysler Improved Numerical Differencing Analyzer 
for Third Generation Computers) was chosen for the basic thermal analysis, and several 
subroutines were written to tailor it for the needs of the particular problem. A portion 
of the tailoring included a program called the pre-preprocessor, which generated the 
data required as input for the CINDA-3G program. 

This pre-preprocessor and allied subroutines comprise a package to obtain a para- 
metric thermal analysis of any multilayer covering on a spheroid. As a result, a mini- 
mum amount of input is required to specify the model geometry so that perturbations in 
the geometry configuration may be introduced with a minimum of effort. 


INTRODUCTION 

The long-term storage of liquid hydrogen in space can be effectively accomplished 
by using high-performance multilayer Insulation (MLI) systems. A typical system would 
use a large number of highly reflective shields (aluminized polyester film) separated by 
low-thermal-conductivity spacers (e. g. , silk netting). 

A cross-sectional view of a sun-oriented vehicle with an unshrouded liquid-hydrogen 
tank covered with MLI is shown in figure 1. The payload intercepts all solar energy and 
thus represents the only high-temperature heat source in the vehicle. The direction of 
the anticipated heat flows could be as indicated by the arrows in figure 1. The only heat 



being transferred into any of the tanks is tlirough the MLI on the top half of the liquid- 
hydrogen tank. Therefore, this discussion will consider the thermal performance of the 
MLI on the top half of the liquid- hydrogen tank only. 

This analysis will be general for a steady- state study of any oblate spheroidal tank 
and payload geometry. However, to simplify the discussion, the programs herein will 
be described with respect to a particular problem; that is, a steady-state analysis of a 
variable-thickness shield of MLI on the top half of an oblate spheroidal tank viewing a 
constant-temperature payload. A complete discussion and the results of this problem 
for the geometries considered may be found in reference 1. 

A cross-sectional view of the variable- thickness MLI system on the top half of the 
liquid- hydrogen tank is shown in figure 2. The net heat flow by thermal radiation into the 
top surface of the MLI is transferred within the MLI by 

(1) Conduction in the lateral direction parallel to the shields, 

(2) Combination of radiation and conduction in the normal direction perpendicular to 

the shields, 

For a steady-state condition to exist, the net thermal radiation received by the top 
surface of the MLI must be exactly equal to the total heat transferred through the MLI 
and into the liquid- hydrogen tank. The problem then becomes one of solving for the net 
heat flows after temperatures resulting from local heat balances throughout the MLI 
have been established. 

There are a number of computer programs available which, when fed the proper in- 
put data to describe the model of a physical problem, will perform a thermal analysis 
and output the temperature distribution. When these programs are used, both building 
the mathematical model which describes the physical geometry and the preparation of the 
necessary input data are most tedious. To alleviate this difficulty, the computer sub- 
routines described in this report were written to provide a ready means to generate this 
input data. 

These subroutines are specifically designed to solve the problem described in refer- 
ence 1, but can be generally applied to any tank in the shape of an oblate spheroid with 
layers of material covering the tank, or for a shell of varying thickness having the same 
general shape. These computer programs, then, generate the data required by 
CINDA-3G (ref. 2), which is the particular thermal analysis program chosen by the 
authors because of its versatility. The capability to readily effect small changes in 
geometry added to the versatility and afforded a substantial saving in time and effort. 

The routines to be described herein are (1) a pre-preprocessor which takes the input 
data and converts it to the CINDA-3G data as prescribed by that program, (2) a view fac- 
tor program which computes the geometric view factors for radiation in an enclosure, 

(3) an output routine which graphically displays the temperatures, thus making it easier 
to get an overall view of the temperature distribution, and (4) various other subroutines 
that are used for special operations by routines (1), (2), and (3). 
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ANALYSIS 


Thermal Mode! 

The thermal model for this problem is shown in figures 3 and 4. It consists of a 45° 
wedge of the MLI on top of an oblate spheroidal tank. This wedge is divided into a num- 
ber of internal nodes. 

For convenience, the geometry of all the nodes in a given column is assumed to be 
identical in the normal direction, having numerical values for the node width, node 
length, the lateral cross-sectional area (ALAT), and the normal area (ANORM) based on 
the dimensions on the top surface of the MLI. This should introduce no significant error 
when the maximum MLI thickness is small relative to the minor tank radius. 

The node length for each segment (column) is described by specifying the angles 

02» • • • j integers). The depth of all nodes is constant and is 

specified by two input variables: the number of shields in each node, NSPL, and the 
thickness of each shield, DELS. The product of NSPL and DELS then is the depth of 
each internal node. The number of nodes in each column is an input variable, NL, 
which is a function of the MLI thickness being simulated in that particular column. The 
program is set up to accept a maximum of 30 segments and up to 50 nodes in each col- 
umn. For the particular problem under discussion, 10 segments were used with varying 
numbers of layers (see fig. 3). The computed performance of this wedge is then multi- 
plied by 8 to obtain the results for the entire top half of the tank. 

A constant- value boundary node represents the temperature of the tank and the tank 
surface is held at this constant temperature. The payload is simulated by a flat disk 
divided into concentric rings (maximum, 10 rings). Space is simulated by a constant- 
temperature node with temperature at absolute zero and an emissivity equal to imity 
(blackbody). 


Thermal Radiation 

An enclosure composed of the payload, the surface of the MLI, and an imaginary 
surface representing space is used to compute the net thermal radiation received by the 
top surface of the MLI. The enclosure formed with the various surfaces labeled is shown 
in figure 5. The enclosure used in this discussion is composed of the 10 tank surface 
areas, five payload surface areas, and an imaginary surface simulating space, for a 
total of 16 surfaces. 

The payload surface represented by a flat disk is divided into concentric rings of 
approximately equal area (labeled A^j^ to A^g in fig. 5). Each of these rings is consid- 
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ered to be a separate surface of the enclosure, and as such, each ring may have a dif- 
ferent surface emissivity and a different constant temperature which is considered to be 
an infinite heat source. This is then assumed to be a multiple gray surface enclosvire 
and simulates a radiosity network. A CINDA-3G subroutine, IRRADl, is then used to 
compute the net radiant heat flow rates to each surface of the enclosure. Subroutine 
IRRADl uses a method described in reference 3 to obtain the net radiant heat flow to 
each surface. 


Computer Program - General 

CINDA-3G is a two-pass or double-phase operation consisting of a preprocessor 
phase and an execution phase. The user input data (description of the model) to 
C1NDA-3G are converted into FORTRAN IV language subroutines by the preprocessor. 
These subroutines are then passed onto the system FORTRAN compiler for compilation 
and execution for solution of the problem. 

Since the user input data as described in reference 2 are for a fixed geometric 
model, much extra labor would be necessary to change the input for a slightly modified 
geometry. Thus, a set of subroutines were written to accept a general description of 
the model and to generate the user input data so that the geometry can easily be varied 
by small changes in input data. This set of author-written subroutines was added to the 
CINDA-3G system as an additional phase of operation, making a three-pass or three- 
phase operation. This additional phase, called the pre-preprocessor, is executed first 
to generate all the necessary input data as required by the CINDA-3G preprocessor using 
a minimum amount of punched card input data. The CINDA-3G data are described in ref- 
erence 2 (pp. 4. 1 to 4. 21). Control is then passed to the CINDA-3G preprocessor 
(phase 2) to process this generated data and to generate the subroutines for the third 
phase of the operation to perform the execution of the problem. 

Because of this method of combining the three phases of operation and the use of 
certain portions of the computer operating system, this overall program as described is 
usable only on IBM 7094 11/7044 Direct Couple Systems running under the operating sys- 
tem IBSYS version 13, with a FORTRAN IV compiler. However, the CINDA-3G program 
is available for other computer configurations, for example UNIVAC 1108; and the rou- 
tines written by the authors as described in this report could easily be modified to oper- 
ate on these other computers since they are entirely written in FORTRAN IV. 


Finite Difference Equations 

The input to CINDA-3G is the description of the model. This description is given by 
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designating each incremental volume, called a node, by a number and then describing 
the heat paths between nodes by specifying the node numbers of the two adjoining nodes 
and a value which represents the conductance between these two nodes. Conductance is 
the area of heat flow times the thermal conductivity divided by the conductor length. 

Heat flow to boundaries is specified by supplying heat -transfer coefficients or radiation 
constants as the conductance between the two nodes. 

When a network (model) to the problem has been described, the solution is obtained 
by solving a set of finite difference equations, which are an approximation to the partial 
differential equations of the diffusion type 

— =Q!V^T + S (1) 

at 


where T is temperature, t is time, a is dlffusivity 


a 



( 2 ) 


is the Laplacian operator in x, y, z 


coordinates, and S 


is a source term defined 


as 


S = 


qo! 

k 


(3) 


where q is heat rate per unit volume and k is thermal conductivity. 

Since we are interested in only the steady-state solution for this particular problem, 
the time derivative is equal to zero, and equation (1) reduces to Poisson’s equation: 

aV^T + S = 0 (4) 

Substituting equations (2) and (3) into equation (4) results in 


kV^T , qk _ Q 


(5) 


Then since pCp may be assumed to be a constant for steady- state computation, multiply 
both sides of equation (5) by pC to get 
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+ q = 0 


( 6 ) 


Then 


kV'^T + q = k 


a^T ^ a'^T . a^^T 
+ + 


vax 


ay 


az 


kV^T + q ^ 

Ax yax" 


ax\ ^ / 3T 

ax'*‘y Ay \ay“ 


aT \ 

+ — 1 

( aT 

9y7 

+ 1 

Az 

\az" 



( 7 ) 


where the plus and minus signs in the first partial term indicate that they are taken on 
the negative or positive side, respectively, of the point under consideration. If three 
consecutive points (1, 0, and 2) ascending in the x direction are considered, the finite 
difference of the x portion of equation (7) is 


Ax \ ax" 




( 8 ) 


applying this to equation (7) for x, y, and z for the points 

1, 0, 2 in X direction 
3, 0, 4 in y direction 
5, 0, 6 in z direction 

yields: 




( 9 ) 


6 



I 


taking common denominator of volume V = (Ax)(Ay)(Az) and = (Ay)(Az), 

A = (Ax)(Az), a = (Ax)(Ay) and combining equation (9) with equation (6) to obtain 

y 2 




kA 

+ ^(Te-T(,)+i = o (10) 


now 


let Gj = kA^/Ax", Gg = kA^^/Ax , Gg = kAy/Ay“, . . . , etc. Then 

°1 (Ti - To) + ^2 (t 2 - To) °3 (T3 ' Tq) - (T 4 ‘ T„) + G5 (Ts ‘ Tq) 

+ Gg (t^ - T 


g-xg)+q = 0 (11) 




= To(Gi + G2 °3 * °4 * °5 * S) 


or 


_ ^l"^! ^2^2 ^ ^3^3 ■'■ *^4^4 ^ 6^6 ^ 

Tg — ^ - 


Gj + Gg + G3 + G4 + Gg + Gg 


(13) 


or, in more general terms, 


To = 


. B^iTt 


“O 




where i ranges over all neighboring nodes to Tg. 


(14) 
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Since all the temperatures T. in equation (14) are not known, assume some initial 
temperatures and solve this equation for Tg; when this is done over all the nodes in the 
network, a new temperature distribution will be obtained. Iterate in this fashion with 
equation (14) over all the nodes in the network until the temperatures obtained do not 
change between two consecutive iterations by more than some prescribed value, called 
the convergence criterion. This is the method used by the execution routine CINDSS (see 
ref. 2). 


Geometric Radiation View Factors Between Surfaces on a 
Flat-Surfaced Payload and a Spheroidal Tank 

The geometric radiation view factor between any two surfaces and A. shown in 
the following sketch is defined by equation (15) 



(b) 


«/a. 


cos cos <p^ dA. dAj 
ttR^ 


(15) 


8 



where F. . represents the portion of the total energy, either emitted or reflected by 
J 

surface i, that is intercepted by surface j. As shown in the sketch, and are the 

angles between the respective normals to the surface elements dA^ and dA^ and the 
connecting line R. 

The double integral in equation (15) is approximated by the following double summa- 
tion in order to obtain a numerical solution on the digital computer; 




(16) 


The radiosity subroutine used to calculate the net heat into each surface of the enclosure 
requires the following numerical input: 

A,,t area of each surface in the enclosure 

% 

T^^ temperature of each surface in the enclosure 

emissivity of each surface in the enclosure 
E 

AxF area times geometric view factor for each interchange in the enclosure (see 
ref. 2, p. 6. 6, 2) 

The procedure described in this section first determines the individual values of AxF 
for each of the tank surfaces viewing each of the payload surfaces. Once these primary 
values are determined, all other values of AxF are readily determined. The equations 
used to define the other values of AxF are discussed at the very end of this section. 

The view factors, and subsequently the net q into each of the surfaces in the radios- 
ity enclosure, are computed for the entire top half of the tank and the entire payload. As 
described earlier, in the discussion of the model, the top node of each column of MLI is 
a 45*^ segment. These segments are extended completely around the tank to form a ring 
and the net q is computed for each of these rings. These q's are then divided by 8 to 
obtain the net q into each node on top of the MLI. The method of computing the view 
factors described herein further breaks up each of these rings on the MLI and also breaks 
up the payload surfaces, described on page 3, into smaller rings. Then these small 
areas with their respective view factors are added together to obtain the AxF for each 
payload surface and each MLI surface, where each MLI surface is a ring as described 
previously. 
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The actual view factors are calculated as follows: The payload is arbitrarily divided 
into concentric rings each having a width of 0. 1 foot (3. 05 cm). Each ring with area A 
is then further divided into incremental areas aA^. As shown in figure 6, the subscript 
p will pertain to the payload surface, and the subscript t will pertain to the tank sur- 
face. The number of incremental areas in each ring Nj is varied so that the numerical 
value of AA remains approximately constant over the entire payload. Thus, for the 

Jr 

payload ring shown in figure 6, Nj is the nearest whole number determined by the fol- 
lowing expression: 


Nj = 207TRp (17) 

The number of these rings of 0. 1-foot (3. 05-cm) width which are to be included in each 
of the individual payload surfaces is then readily determined. 

The MLI surface on the tank is arbitrarily divided into 90 increments of 1° each, 
measured from the vertical z-axis. As shown in figure 6, the incremental area extends 
in a ring completely around the tank. The location of this ring is specified by the 

angular distance from the vertical z-axis. Each ring with area A^ is further divided 
into incremental areas aA^ of constant size over the entire tank surface. The number 
of individual rings to be included in each of the separate tank surfaces used in the en- 
closure is then readily determined once the angular size of each of the surfaces is speci- 
fied. 

The variable MLI thickness over any particular tank surface area is accounted for by 
reading in the number of nodes used at that particular location. The numerical values of 
a and b (the major and minor semi-axes of the spheroid formed by the surface of the 
MLI) are determined individually for each of the tank areas used. 

a ~ + AS 

b = + AS 

where AS is the thickness of the MLI and and Y„ are the uninsulated tank dimen- 

r r 

sions. Thus for any particular surface area of the tank, the existing MLI thickness over 
that particular area is considered to exist over the entire tank surface. 

Both the tank and the payload are symmetrical about their respective z-axes. Be- 
cause of this symmetry and the fact that the z-axes of both coordinate systems coincide, 
the geometric view factor of any one element aA^. along a given tank ring to any given 
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complete payload ring is independent of position. Thus, equation (16) may be written in 
the following manner: 


Vt,p 




(19) 


where and Ap are the tank and payload ring areas, respectively. Thus, the view 
factor F. of any ring on the tank to any ring on the payload is the summation of the 
view factors to all the incremental areas making up the payload ring. 

To simplify the geometry, the typical element AA^ on any given tank ring, selected 
to view all elements AA on any given payload ring, is located in the yz plane. In this 
plane, the equation of the MLI surface is 



( 20 ) 


The following geometric relations are then valid; ' 


Rj. = ab 






2 . 2 , 

sin xj/ 


m 


9 9 

a^ cos 


m 


yt = ^t 

Zj = cos V'jn 
Ah = b - z^ 


( 21 ) 


For each value 
areas AAp. Thus, 
will be a function of 
ure 6: 


of used, the entire payload will be divided into incremental 

Rp will be a function of the payload ring being considered, and 

the location of AA . The following relations can be seen from fig- 
P 
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where Nj is determined by equation (17). Differentiating equation (20) gives the slope 
at any point on the curve: 


dz _ y 
dy ^2 z 


(23) 




Since the payload surface is in the xy plane, 




(26) 


(27) 
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The vector R. from aA. to AA can be expressed as 
t,p t p 

^t,p = (""p - + % - yt> i + (^p - ^t) ^ 


( 28 ) 


Since x^. = 0, 


\p = V'"(yp-yt)^ 


(29) 


R = 




(30) 


Thus the unit vector R. becomes 

t, P 


- X i + (yp - y^) j + (Zp - zp k 

\ p 


R 


and 


(31) 


^p, t - p 


(32) 


The angles cp. and <p 

^ Ir 


can now be determined by forming scalar products: 



(33) 


cos. cp^ = N„ 
P 


R = 

^p,t — 


The tank surface area A of any ring is determined with the equation obtained from the 
following derivation: 
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Thus 


A = 277a 
P 




2 2 1 
w z + — 

2w 



(34) 


Therefore all quantities required in equation (19) are determined. 

To obtain the value of between specific areas m and n is simply an 

appropriate summation of the results of equation (19): 


A F 
m m, n 


E 

All tank 
rings in the 
m“^ segment 


E '"‘'■t.p 

All payload 

rings in the 
th 

n ^ segment 


(35) 
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All other values of can now be generated using the following relations: 

(1) Since no surface "sees" itself, m “ ® 

(2) Since no two surfaces on the tank "see" each other, F^ n ~ sur- 

faces m and n. 

(3) Since no two surfaces on the payload "see" each other, F^^ n ” ® P^y~ 

load surfaces m and n. 


All surfaces 

(4) By definition, F 


m, n 


= 1 . 


(5) By reciprocity, A^F^^ = A^^F^^ 


DESCRIPTION OF SUBROUTINES 

As previously discussed, the program is divided into three phases: (1) the pre- 
preprocessor, which by using the input data described in appendix A of this report com- 
putes the input data needed to run a problem with CINDA-3G and generates a tape accord- 
ing to the formats given on pages 4. 1 to 4. 2 of reference 2; (2) the preprocessor, as de- 
scribed in reference 2; and (3) the execution phase, which actually computes the temper- 
ature distribution. 


Pre-preprocessor Phase 

The pre-preprocessor is composed of four subroutines which read data cards de- 
scribing the geometry, compute the data required by the CINDA-3G system, and write 
this data onto a tape to be used in phase 2. 

The subroutines in the pre-preprocessor are as follows: 


Subroutine 

name 

FORTRAN 
call name 

Brief description 

PREPRE 

Main program 

Reads the input data, computes the data required 
by CINDA-3G 

VIEWFC 

VIEWF 

Computes the geometric view factors for the ra- 
diation interchange used in subroutine IRRADl 

RAYOUT 

RITARY 

Called by PREPRE; writes arrays of numbers 
onto the input tape as prescribed by CINDA-3G, 
p. 4. 10 of ref. 2 

COMDAT 

BLOCK COMMON 

Block common subroutine 
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Preprocessor Phase 


The preprocessor is the second phase of the operation. A list of the subroutines and 
a description of each is given in reference 2 (pp. 8. 1 to 8. 3). The only variation is that 
the preprocessor reads the data from the tape that was prepared by the pre-preprocessor 
phase instead of from cards. 


Execution Phase 

The actual computing to determine the temperature distribution is done in the execu- 
tion phase. In this execution phase the following subroutines are used: LINKO, EXECDK, 
VARIDK, VAR2DK, OUTDK, which are generated by the preprocessor phase of 
CINDA-3G. 

LINKO is the main program. Its only fimction is to call the INPUT routines to read 
in the CINDA input data and call EXECDK. 

EXECDK is called the execution block of the program. It retains control until the 
problem is solved. Subroutines VARIDK, VAR2DK, and OUTDK are called by EXECDK 
or by other subroutines called from EXECDK, 

VARIDK, VAR2DK, and OUTDK are the subroutines generated from the VARIABLES 
1 block, VARIABLES 2 block, and OUTPUT block, respectively, that are described in 
reference 2, The order of calculation in the execution phase of the program is described, 
along with a flow chart, in reference 2 (pp, 3, 12 and 3, 13). 

Several subroutines were written by the authors and incorporated into the execution 
phase in the user subroutine section. These subroutines may be easily modified or re- 
written by the user because they must be included into the program deck as FORTRAN 
source cards. The subroutines are 


Subroutine 

name 

CNDCTR 

OUTEND 

DISTRT 

GRAFIC 

PUNTEM 


FORTRAN 
call name 

GETG 

ENDOUT 

TDIST 

PICTR 

PCHT 


Subroutine CNDCTR computes the value of all internal conduction connectors G^ 
between two nodes i and j. 
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( 36 ) 



A 


k 




J 


where 

A area of heat flow 
k thermal conductivity 

I distance from the center of node i to center of node 3 
If the thermal conductivity varies with temperature 

a = (37) 

k, kj 

where k^^ and k^ are temperature -varying thermal conductivities at nodes i and j. 
These conductors are computed just prior to each iteration of the network solution. 

Subroutine CNDCTR as listed in this report was written to permit the use of very 
high thermal conductivity shields (pure aluminum as opposed to aluminized mylar) to in- 
crease the lateral thermal conductivity of some layers of the MLI, The nodes in these 
layers with high conductivity will be denoted as hi-k nodes. All other interior nodes with 
a lower conductivity will be denoted as standard nodes. It should be noted that this is 
conduction in the lateral direction only. The thermal conductivity for the standard nodes 
is determined by straight-line curve fit. Input variables CONI and,CON2 are the a 
and b, respectively, for the line k = aT + b. The data for the hi-k nodes are input in a 
table with variable name AKVST. The format for this table is found in appendix E. For 
the data used in the sample problem for both the standard nodes and the hi-k nodes see 
figure 8 of reference 1. The thermal conductivity for conduction in the normal direction 
^effective ^ combination of conduction and radiation (see ref. 1 for details). Sub- 
routine CNDCTR also calls subroutine IRRADl to obtain the net Q into each surface of 
the radiosity enclosure. 

Subroutine OUTEND calls subroutine GRAFIC and then computes and prints out the 
heating rates Q of the fuel tank to MLI, QNORM, and also the net heat into the MLI on 
the top surface of the MLI, QTS. QNORM and QTS are printed in tabular form to show 
the Q in or out of each segment on the tank. Further, these values are then multiplied 
by 8 to obtain the Q in and out for the entire tank. The names of these output variables 
are QN 8 and QTS 8 . These printouts occur each time the print switch (NPRINT) calls for 
an output of the temperature distribution. 
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Subroutine TDIST reads the initial temperature distribution from tape unit SYSUT5 
which was written in the pre -preprocessor. 

Subroutine GRAFIC is a special subroutine which prints out the temperature distri- 
bution pictorially. Blocks are made, simulating the nodes, and are printed to simulate 
the shape of the shield. Inside each block is printed a node number and the temperature 
of that node. A node number with a minus sign in front denotes this node is a hi-k node, 
as described in the section on subroutine CNDCTR. If the top of the box is made up of 
dollar signs ($$$$$$) instead of asterisks (**+***) it denotes a surface with a high sur- 
face emissivity (see sample output listing on p. 31). 

Subroutine PUNTEM punches cards with the current temperatures if either the num- 
ber of iteration loops or the amount of computer time requested on the TCP card is ex- 
ceeded. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 2, 1970, 

129-04. 
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APPENDIX A 


COMPUTER PROGRAM SYMBOLS 


FORTEIAN 

name 

Corresponding 

mathematical 

symbol 

Units in 
program 

Description 

A 

A 

norm 

ft2 

Conduction area in normal di- 
rection 

A FA 


ft2 

th 

Product of i*"" surface area 
and the geometric view factor 

XU AU 

between i“ and j surfaces 

AKVST 

k 

Btu/(hr)(ft)(°F) 

Thermal conductivity table 

ARCLN 

— 

ft 

Arc length 


AREAL 


ft2 

Conduction area in lateral di- 
rection 

BOTLIM 

— 

— 

Integrand of equation (34) eval- 
uated at lower limit 

BTA 


deg 

Angle measured from minor 
axis- to segment boundary on 
tank 

C 

— 

— 

Array of user constants to 
CINDA-3G 

CINCON 

— 

— 

Array of control constants to 
CINDA-3G 

CONNAM 

— 

— 

Array of control constant 
names 

CONI 

~ — 

“ “ 

Constant = a^ 

a straight-line 
curve fit for k 
for lateral con- 

CON2 

— 

— 

Constant = b^ 

duction; k = aT-i-b 

DAMPA 

— 

— 

Damping factor (CINDA control 
constant) 

DAREA 

At 

ft2 

Incremental area on MLI for 
view factor computation 
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FORTRAN 

name 

Corresponding 

mathematical 

symbol 

Units in 
program 

Description 

DELDMP 

— 

— 

Delta to add to change damping 
factor 

DELS 

— 

ft 

Thickness of each shield 

DMLI 

AS 

ft 

Thickness of MLI at a given 
segment 

DPTH 

— 

ft 

Thickness (depth) of each node 

EMISP 

€ 

— 

Emlsslvlty of each surface of 
payload 

EMISS 

e 

— 

Emlsslvlty of space 

EMIST 

€ 

— 

Emlsslvlty of each surface on 
MLI 

G 

G 

Btu/(hr)(°F) 

Conductor value for heat trans- 
fer 

H 

h 

ft 

Tank-to-payload spacing 

lEXECN 

— 

— 

Name of CINDA execution 
routine to be used 

LCK2 

— 

— 

Tape unit SYSCK2 CINDA Input 
written on 

LUTS 



Tape unit SYSUT3 used as tem- 
porary storage for FORTRAN 
subroutines to be used In exe- 
cution phase 

LUTS 



Tape unit SYSUT5 Initial tem- 
peratures stored on this unit, 
to be read In execution phase 

MAXNL 

— 

— 

Maximum number of nodes at 
any segment 

NDAMP 

— 

— 

Change damping factor every 
NDAMP Iterations 

NHIK 

— 

— 

Number of hlgh-conductlvlty 
layers 
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FORTRAN 

name 

Corresponding 

mathematical 

symbol 

Units in 
program 

Description 

NIN 

— 

— 

Tape unit, input unit 

NL 



Number of nodes in MLI in 
normal direction for each seg- 
ment 

NNODES 

— 

— 

Total number of nodes in con- 
figuration 

NNPAY 

— 

— 

Number of nodes on payload 

NNSEG 

— 

— 

Number of segments on tank 

NOFA 

— 

— 

Number of entries in view fac- 
tor matrix 

NOUT 

— 

— 

Tape unit, output unit 

NPAY 

— 

— 

Same as NNPAY above 

NPRINT 

— 

— 

Print temperatures every 
NPRINT iterations 

NS EG 

— 

— 

Same as NNSEG above 

NSPL 

— 

— 

Number of shields in each node 

NSTMTS 

— 

— 

Number of FORTRAN cards 
read in as user program to go 


into SUBROUTINES block of 
CINDA-3G 


NSURF 



Number of surfaces in radiation 
enclosure 

NTAB 

— 

— 

Total number of entries in 
K-vs-T table (AKVST) 

NUMHIK 



Array containing layer num- 
bers of high- conductivity lay- 
ers 

PSI 

\L 

^m 

— 

Angle from top of tank to Rj^ 

QNORM 

Q 

Btu/hr 

Heat flow from MLI into tank 
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FORTRAN 

Corresponding 

Units in 

Description 

name 

mathematical 

symbol 

program 


QTS 

Q 

Btu/hr 

Heat flow from enclosure into 
top of MLI 

R1 


ft 

Distance from tank center to 
incremental surface i 

RIX 


ft 

Distance from tank center to 
incremental surface i + 1 

RP 


ft 

Radius of payload 

SLAY 



Number of nodes in MLI in 
normal direction at each 1^ 
increment 

SPACES 

e 

— 

Emissivity of space 

TCARDS 



Switch: If equal to zero, use 
temperatures input with input 
data; then distribute initial 
temperatures throughout con- 
figuration with subroutine 
TDIST. If not equal to zero, 
read initial temperature distri- 
bution from cards. 

THETA 

e 

deg 

Array of angles specifying 
angular segments on tank 

THIKN 

— 

— 

Same as DPTH above 

TPAY 

T 

°R 

Array temperatures of payload 
nodes 

TSPACE 

T 

°R 

Temperature of space node 

TSURF 

T 

°R 

Initial temperature on MLI 
surface 

TWALL 

T 

°R 

Temperature of tank wall 

UPLIM 

— 

— 

Integrand of equation (34) 
evaluated at upper limit 

X 

^r 

ft 

Major semi-axis of tank 
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FORTRAN 

Corresponding 

Units in 

Description 

name 

mathematical 

symbol 

program 


XMl 


— 

Slope 

Y 

Yr 

ft 

Minor semi-axis of tank 
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APPENDIX B 


FLOW DIAGRAM 


General program 


Pre-preprocessor 


START START 
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APPENDIX C 


DECK SETUP 

The following is the deck setup to run this program at Lewis Research Center on the 
IBM 7094-7044 Direct Couple System using IBSYS version 13: 


ccl 

cc8 

ccl6 

$ID 

JO 

USER NAME 

$TCP 


TIME=XX, PAGES=YY 

$SETUP 

10 

CIN3G 


$IBJOB 

FORTRAN IV source or binary decks of pre-preprocessor as described in this report. 
(Includes subroutines PREPRE, VIEWFC, RAYOUT, and COMBAT.) 


$DATA 

FORTRAN IV source cards of user-written subroutines. For this problem these sub- 
routines include CNDCTR, OUTEND, DISTRT, GRAFIC, and PUNTEM. Any other 
user -written subroutines may be added here. 

ENDX 

Input data as described in appendix 

$IBSYS 
$REWIND 
$SWITCH 


SYSCK2 

SYSINl, SYSCK2 
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APPENDIX D 


TAPE USAGE 

Tape units must be made compatible between CINDA-3G and the computer system. 
A table of tape usage presently being used at Lewis Research Center only with the 
CINDA-3G program is included to aid in making this compatible with other systems. 


SYSUNI 

name 

FORTRAN 
number in 
CINDA-3G 

Variable 
name in 
program 

Lewis 

system 

number 

Function 

UT3 

2 

LUT3 

14 

Copy of original prob- 
lem data 

UT4 

3 

LUT4 

15 

Parameter change data* 

UT8 

4 

LUTl 

3 

Data number definitions 

INI 

5 

NIN 

5 

Input 

OUl 

6 

NOUT 

6 

Output 

UT2 

8 

LUT2 

13 

NA-NB pairs, data 
number definitions for 
parameter changes 

CKl 

9 


10 

CINDA-3G master tape 

CK2 

10 


16 

New master tape if up- 
dating, also used as 
problem recall data 
tape* 

UT5 

11 


7 

Problem store data 
tape* 

UT6 

12 

LBD3 

1 

Data tape (original 
problem and all param- 
eter changes) 

PP1(LB4) 

13 


Calcomp 

System plot tape 

LB3 

14 

LB4P 

9 

Program tape (contains 
generated FORTRAN 
routines; LINKO, 
EXECTN, VARBLl, 
VARBL2, OUTCAL) 
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SYSUNI 

FORTRAN 

Variable 

Lewis 

Function 

name 

number in 
CINDA-3G 

name in 
program 

system 

number 


UT7 

15 

LUT7 

2 

Variables 1 calls gen- 
erated from node and 
conductor data blocks* 

UT9 

16 


4 1 

I 

r Internal file for reread 

UT9 

99 


^ J 

I 

LB2 

17 


8 

Overlay tape 


The tapes marked by an asterisk need not be assigned if the particular options are not 
used. The STOREP option in CINDA-3G requires assigning and saving tapes on units LB3 
and UTS, The RECALL option requires assigning and mounting these two tapes onto 
units LBS and CK2, respectively. 
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APPENDIX E 


DESCRIPTION AND FORMAT OF INPUT DATA TO PRE-PREPROCESSOR 

Input Data Description 


CARD 1 


FORMAT( 

NSEG 

NPAY 

NS PL 

DELS 

TWALL 

TSURF 

TSPACE 

EMISS 

SPACEE 

RP 


14, 213, 7F10.5) 

NO. OF SEGMENTS (MAXIMUM 30) 

NO. OF SURFACES, (SEGMENTS) ON PAYLOAD (MAXIMUM 10) 
NO. OF SHIELDS PER LAYER 

THICKNESS OF EACH ( SHI ELD^-SP ACER ) FT. 

TEMPERATURE OF TANKWALL DEG.R 
TEMPERATURE OF TOP SURFACE OF SHIELD DEG.R 
TEMPERATURE OF SPACE DEG.R 

EMISIVITY USED IN EQUATION FOR K(EFF), 

CONDUCTION IN NORMAL DIRECTION 
EMISIVITY OF SPACE 
RADIUS OF PAYLOAD FT. 


CARD 2 FORMAT! 

TCARDS 


X 

Y 

H 

CONI 

COM2 

lEXECN 

CARD 3 FORMAT! 
NPRI NT 
NDAMP 
OELDMP 


6F10.5» A6) 

SWITCH- TO READ INITIAL TEMPERATURES 
IF TCARDS =0, CARDS WILL NOT BE READ. THE PROGRAM 
USES TWALL AND TSURF AND DISTRIBUTES THE TEMPERATURES 
BETWEEN THESE TWO VALUES AS THE INITIAL TEMPERATURE 
DISTRIBUTION. 

IF TCARDS .GT. 0 THE INITIAL TEMPERATURE DISTRIBUTION 
IS READ FROM CAROS, (DEG. R>. SEE CARDS 11 BELOW. 
LENGTH OF SEMI-AXIS OF ELIPSE, ( TANKWALL ), ALONG X AXIS 
LENGTH OF SEMI-AXIS OF ELIPSE, ( TANKWALL ), ALONG Y AXIS 
TANK TO PAYLOAD SPACING FT. 

CONSTANTS USED FOR STRAIGHT LINE CURVE FIT TO COMPUTE 
THERMAL CONDUCTIVITY FOR G VALUES FOR LATERAL 
CONDUCTION OF STANDARD NODES. K=C0N1*T*C0N2. 

NAME OF EXECUTION ROUTINE CINDSS OR CINDSL 

2110, F10.5) 

DELTA PRINT, PRINT TEMP. DIST. EVERY NPRINT ITERATIONS 
CHANG DAMPING FACTOR EVERY NDAMP ITERATIONS 
DELTA TO ADD TO DAMPING FACTOR 


CARO 4 FORMAT! 18, 9F8. 3) 

CINCONU) CINDA CONTROL CONSTANTS, AS REQUIRED BY EXECUTION 
ROUTINE BEING USED. CONSTANTS INPUT IN ORDER AS 
FOLLOWS ACCORDING TO FORMAT GIVEN ABOVE, 

NLOOP, ARLXCA, OAMPA, ORLXCA, OAMPO, TIMEO, 

TIMENO, OUTPUT, OTIMEI, CSGFAC. 

CARDS 5 NTAB, ! AKVST ! I ) , I= 1 ,NTAB J FORMAT ( 16 , 12F6. 1/ ! 6X , 12F6 . 1 ) ) 

TABLE OF TEMPERATURE VS. THERMAL CONDUCTIVITY FOR 
LATERAL CONDUCTION OF HIGH CONDUCTIVITY LAYERS. 

NTAB= TOTAL NO. OF ENTRIES IN TABLE (TEMPS AND K'S) 

data form NTAB TEMPI K1 TEMP2 K2 TEMPN KN 

UNITS — K BTU/HR-FT-DEG.R 
T DEG.R 

CARDS 6 NL NUMBER OF LAYERS IN EACH SEGMENT FORMAT ! NI2) 



CARDS 7 THETAO ANGLE OF EACH SEGMENT* IN DEGREES* IF ALL SEGMENTS 

ARE EQUAL INPUT ONLY THE FIRST. FORMAT (N 12) 
THE SUM OF THE N ANGLES MUST EQUAL 90 

CARDS 8 ON TWO CAROS (10F8.5) 

EMISIP(I) EMISIVITY OF EACH PAYLOAD SEGMENT f MAXIMUM 10) 
TPAYU) TEMPERATURE OF EACH PAYLOAD SEGMENT IMAXIMUM 10) 

CARDS 9 EMIST EMISSIVITY OF EACH SURFACE* RADIATION (8F10.5) 

CARD 10 NHIK NO. OF HIGH CONDUCTIVITY LAYERS (3012) 

NUMHIK LAYER NUMBERS WHICH ARE HI COND. 

CARDS 11 NNODES IF INITIAL TEMPS READ FROM CARDS - NO. OF NODES 
T(I) INITIAL TEMPERATURES* DEG. R ( I 8* 9F 8. 1/10F8. 1 ) 


Format of Input Data 
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APPENDIX F 


INPUT DATA FOR SAMPLE PROBLEM 


SHILL 


IrtG STUDY CDMFIGURATiaN 

m. 24 






10 


5 5 

. 0016667 

37.0 

295.0 

0.0 

0, 

24 ■ - 1 . 0 


5.0 




4.571 

3 

.233 

0.8175 

000047 

.0045 CINDSS 





125 

2000 

0.1 




■ - 



s 

s 

99 , 

0001 0 

.1 








It 


37. 

.8 

40 

.873 

45. .913 

50. 

.970 

55. r.OlO 

60. 

1.045 



62. 

1.055 

64 

1.065 

66. 1.0875 

68. 1 

.105 

70. 1.112 

72. 

1.114 



74. 

1,116 

75 

1.117 

76. 1,116 

78. 1 

.114 

80," 1.112 

82. 

1.105 



84. 

1.090 

86 

1.085 

88. 1.072 

90. 1 

.062 

95. 1.040 

100. 

1.012 



105. 

.985 

no 

.955 

115, -925 

120. 

.900 

125." “,880 

130. 

.857 



135. 

.837 

140 

.819 

145. .800 

150. 

.785 

155.- ,770 

160. 

.757 



165. 

-742 

170 

.730 

175. .720 

180. 

.711 

185.“ .700 

190. 

.695 



1 195- 

,687 

200 

.683 

205. .675 

210. 

.672 

215.' .668 

220. 

. 665 



1 225. 

.662 

230 

.660 

235- .658 

240. 

.655 





C 

201717151311 

8 








£715 

3 

9 6 

6 6 4 4 

4 








.11 


.11 

11 

.11 

.11 






5£Q 


0 520 

.0 520.0 

520.0 

520.0 






. 02^ 



.024 

.024 

. 024 

024 

.024 .900 

- 

900 

-9oq 



.900 









4 dlj 

1 

1315 








1 



OUTPUT FOR SAMPLE PROBLEM 


SHIELDING STUDY CONFIGURATION NO. 24 


» 






!■» 

« 

»* 

» 

* PAYLOAD 1 

* NODE NOS 

NODES 

179 

180 

181 

182 

183 

* 

« 

* SPACE 

* NODE NO. 184 

«■ 

* 

• 

* TEMPS= 

520.00 

520.00 

520.00 

520.00 

520.00 

» 

« 

* TEMP= 0. 

« 

• 

* EMISIV= 

0.1100 

»«•««««»» 

0.1100 

0.1100 

0. 1100 

0.1100 

* 
I- • 

* EMISIV 1.0000 

• 

»« 

ANGLE OF EACH 
SEGMENTIDEG. ) 

27.0 

15.0 

9.0 

9.0 

6.0 

6.0 

6.0 

4.0 4.0 

4.0 

EMISIVITY OF 
EACH SEGMENT 

0.02AO 

0.02A0 

0.02A0 

0.0240 

0.0240 

0.0240 

0.9000 

0.9000 0.9000 0 

.9000 


TANK TO PAYLOAD SPACING^ 0.8175 FT. 
LAYER 
NO. 


24 

• 

25 • 



• 

340.27* 


23 

« 

24 * 



• 

331.26* 


22 

• 

23 * 

48 • 


• 

321.44* 

317.63* 

21 

«■ 

22 * 

47 • 


» 

310.72* 

306.77* 


31 



20 


• 21 • 46 • 69 • 90 • 

• 298.88* 294.54* 286.25* 262.06* 


19 • 20 • 45 • 68 • 89 * 

• 285.58* 280.90* 272.16* 252.90* 


18 • 19 • 44 * 67 • 88 • 

• 270.27* 265.30* 256.60* 240.99* 


17 

• 

18 • 

43 • 

»«»****»^ 

66 • 

87 * 

108 • 

126 * 






• 

252.02* 

246. 77» 

238.48* 

225.46* 

215.58* 

195.77* 





16 

• 

M S U wk 9 S W fl ( 

17 • 

42 • 

65 • 

86 * 

107 • 

125 • 






• 

228.92* 

223.24* 

215.80* 

205.86* 

197.55* 

186.02* 








»•«»»«»« 4 








15 

• 

“16 ♦ 

-41 * 

-64 • 

-85 • 

-106 • 

-124 • 

-142 • 





• 

195.88* 

188.96* 

183.10* 

178.40* 

173.98* 

170.70* 

167.44* 




14 

• 

15 * 

40 • 

63 • 

84 * 

105 * 

123 * 

141 • 





* 

176.08* 

170.70* 

166.28* 

162.63* 

159.38* 

156.82* 

154.65* 















13 

• 

-14 • 

-39 • 

-62 • 

-83 • 

-104 • 

-122 • 

-140 • 

-156 • 




• 

1^6.57* 

143.90* 

141.98* 

140.39* 

138.87* 

137.60* 

136.23* 

134.96* 



12 

• 

13 * 

38 • 

61 • 

82 • 

103 • 

121 • 

139 * 

155 • 




• 

131.96* 

130.02* 

128.59* 

127.33* 

126.17* 

125.16* 

124.16* 

123.49* 




*4 










11 

• 

-12 • 

-37 • 

-60 • 

-81 • 

-102 • 

-120 • 

-138 • 

-154 • 

-168 • 



• 

110.36* 

109.60* 

108.99* 

108.42* 

107.84* 

107.29* 

106.67* 

106.05* 

105.40* 


10 

• 

11 • 

36 • 

59 * 

80 • 

101 • 

119 • 

137 • 

153 • 

167 • 



» 

132.06* 

101.46* 

100.94* 

100.45* 

99.96* 

99.51* 

99.04* 

98.65* 

98.51* 


9 


10 * 

35 • 

58 • 

79 • 

100 • 

118 • 

136 • 

152 • 

166 • 



• 

91.32* 

90.88* 

90.47* 

90.09* 

89.72* 

89.38* 

89.04* 

88.79* 

88.72* 













8 

« 

-9 * 

-34 • 

-57 • 

-78 • 

-99 • 

-117 * 

-135 * 

-151 • 

-165 • 

-177 • 


• 

74.74* 

74.55* 

74.37* 

74.18* 

73.98* 

73.76* 

73.51* 

73.20* 

73.00* 

72.81* 

7 

» 

8 * 

33 • 

56 * 

77 • 

98 • 

116 • 

134 • 

150 • 

164 • 

176 ♦ 


» 

72.04* 

71.89* 

71.74* 

71.60* 

71.45* 

71.31* 

71.15* 

71.08* 

71.03* 

70.99* 

6 

«■ 

7 * 

32 * 

55 • 

76 • 

97 • 

115 • 

133 • 

149 • 

163 ♦ 

I V ft ft ft ft ft ft 

175 • 


• 

59. 13* 

69.02* 

68 .91* 

68.80* 

68.69* 

68.65* 

68.60* 

68.57* 

68.54* 

68.53* 

5 


6 • 

31 * 

54 • 

75 • 

96 * 

114 • 

132 • 

148 • 

162 • 

174 • 


«• 

55.93* 

65.85* 

65.78* 

65.71* 

65.65* 

65.63* 

65.61* 

65.60* 

65.59* 

65.58* 

4 

« 

5 • 

30 • 

53 • 

74 • 

95 • 

113 * 

131 * 

147 • 

161 • 

173 • 


• 

52.27* 

62.22* 

62.18* 

62.14* 

62.11* 

62.10* 

62.09* 

62.08* 

62.08* 

62.08* 

3 

• 

4 * 

29 * 

52 • 

73 • 

94 • 

112 • 

130 * 

146 * 

160 • 

172 • 


« 

57.90* 

57.87* 

57.84* 

57.83* 

57.82* 

57.81* 

57.81* 

57.80* 

57.80* 

57.80* 

2 

• 

3 * 

28 • 

51 • 

72 • 

93 • 

111 • 

129 • 

145 • 

159 • 

171 • 


• 

52.30* 

52.28* 

52.27* 

52.26* 

52.26* 

52.26* 

52.25* 

52.25* 

52.25* 

52.25* 

1 

• 

2 • 

27 • 

50 • 

71 • 

92 * 

110 • 

128 • 

144 • 

158 • 

170 * 


» 

43.98* 

43.97* 

43.97* 

43.97* 

43.96* 

43.96* 

43.96* 

43.96* 

43.96* 

43.96* 

TANKWALL 


1 

26 

49 

70 

91 

109 

127 

143 

157 

169 

SURFACE 


37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 


* • 


• TANKWALL * 

• NODE N0.= 178 • 

• TEMP.= 37.00 • 
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SHIELDING STUDY CONFIGURATION NO. 24 


NO. SEG 

QTS 

QNORM 

QTS*8 

ONORM*8 

1 

0.06936889 

0.00018959 

0.55495113 

0.00151672 

2 

0.10318061 

0.00029260 

0.82544488 

3.00234082 

3 

0.09149027 

0.00025254 

0.73192215 

0.00202030 

4 

0.04330576 

0.00032500 

0.34644613 

0.00260004 

5 

0.04622486 

0.00023741 

0.36979891 

0.00189930 

6 

0.00980636 

0.00027727 

0.07845092 

3.00221822 

7 

-0.20052147 

0.00028314 

-1.60417175 

0.00226517 

8 

-0.09447652 

0.00018922 

-0.75581217 

0.00151382 

9 

-0.04718843 

0.00019050 

-0.37750745 

0.00152403 

10 

-0.01882567 

0.00018001 

-0.15060541 

0.00144006 

SUMMATION 

0.23647E-02 

0.24173E-02 

0.18917E-01 

0. 19339E-01 
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APPENDIX G 


FORTRAN LISTING OF AUTHOR-WRIHEN SUBROUTINES 

SIBFTC PREPRE 

COMMOM/BCD/ NBC0(60) 1 

COMMOM /TAPE/NIN,NOUT,LCK2 2 

COMMOM /RADIO/ NSEG, NPAY ,NSURF, X, Y, H,RP , THIKN» PIE t SLAY (90) . 3 

1 THETA0(30),A(41),AFA(1681) 4 

DIMENSION 10(14) ,BTA(31) .ERAOt31),THETA(30)tF(30),lTEST(50)*NL(30) 5 

DIMENSION ARCLN(30)tAREAL(30),NSNOD(30)«NUMHIK(30)tC(36) 6 

DIMENSION AKVST(200), EMIST(41), TPAY(IO). EMISP(IO) 7 

DIMENSION ISTMT(14), T(1545), CONNAM(lO)t CINCON(IO) 8 

DIMENSION IFMTK7), IFMT4(4), IFMTHI5) 9 

DIMENSION Zl(9),Z2(9),Z3(9),Z4(9),Z5(9).Z6(9)tZ7(9),Z8(9) 10 

DIMENSION Z9(9).Z10(9).Z11(9) 11 

EQUIVALENCE (NNSEG,C ) , (NNP AY»C( 2 ) ) » (NSURFS.C ( 3 ) ) . ( NNOOES ,C ( 4 ) ) . 12 

X (MAXNL.CIS) ), (NPRINT,C(6) )» (NDAMP,C(7) ) ,(NNLtC(8) ) , I NAREAL t C ( 9) ) . 13 

X (NAREANfCdO) ), (NARCLN,C( 11 ) 1 1 (NTHETA,C( 12) ) , (NITESTtCI 13) ) , 14 

X (NKHlC,C(14)).(NEMISS«C(15))t(NAFA,C(16)),(NSNT,C(17))t 15 

X (SIGMA, C(21 ))> (DPTH»C(22) ), (DELS. C(23) ).(EMISS,C(24) ) . 16 

X (C0N1,C(25) ),(C0N2,C(26)),(DELDMP,C(27)), 17 

X (TABS,C( 29) ) , (HIKCON,C( 30) ) ,(PSWTCH,C(31) ) , (RAD,C(32) ) 18 

X ,(HH,C(33)) 19 

DATA IFMTl(l)/ 42H( 12/ F2 . 0/ ( 10F8 .5 ) ) / 20 

DATA IFMTH(l)/ 30H ( 7X , 17HBCD 3THERMA ,56X ) / 21 

DATA :0NNAM(1)/60H NLOOPARLXCA DAMPAORLXCA DAMPD T IMEOTI MENDOUTPUT 22 

IDTIMEICSGFAC / 23 

DATA ISL.IWBK.IBACK.IVARB.ISPS.ILPS.IENDX 24 

1 /6HCINDSL,6HCNFWBK,6HCNBACK,6HCNVARB,6HL SPCS,6HL LPCS.6H ENOX / 25 

DATA Ml, ONE, ONEM, TEN20, PI, TABS, SIGMA 26 

1 /I, 1.0, -1.0, 6H1.E+20, 3.14159265, 460.0, 0.173E-8 / 27 

DATA 28 

1 ZKl)/ 54H NL, NO. OF LAYERS IN EACH SEGMENT / 29 

2,Z2(1)/ 54H AREAL, AREA, LATERAL CONDUCTION, EACH SEGMENT (HORIZ)/ 30 

3,23(1)/ 54H A, AREA, NORMAL CONDUCTION, EACH SEGMENT (VERT) / 31 

4.2411)/ 54H ARCLN, ARC LENGTH, LENGTHS FOR LATERAL CONDUCTION / 32 

5, Z5(1)/ 54H THETA (DEG.), ANGLE OF EACH SEGMENT / 33 

6, Z6(1)/ 54H ITEST, NEGATIVE NO. IS HI-COND. LAYER, LATERAL CONDCT/ 34 

7, Z7(1)/ 54H AKVST, K(BTU/HR.FT) VS T(DEG R) FOR HI-K LAYERS / 35 

8, Z8(1)/ 54H EMIST, EMMISIVITY OF SHIELD, RADIATION. / 36 

9, Z9 (D/54H AF, AREA*VIEWFACTOR FOR RADIOSITY, FOR EACH SEGMENT / 37 

A,Z10(1)/54H TABLE OF SURFACE NODE NUMBERS, SURFACE EACH SEGMENT / 38 

PIE=PI 39 

PI04= PI/4. 40 

NIN= 5 41 

N0UT= 6 42 

LUT3=2 43 

LCK2=10 44 

LUT5= 11 45 

REWIND LCK2 46 

REWIND LUT3 47 

REWIND LUTS 48 

DO 12 1=1,10 49 

12 CINCONI I )=0. 50 

C 51 

C WRITE EXECUTE, ID, IBJDB CARDS, THEN ALL CINDA3G PREPROCESSOR 52 
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CARDS OVITO TAPE SYSCK2 (9) 53 

5A 

WRlTE(LCK2t2000) 55 

56 

57 

READ PRE-PREPROCESSOR DATA AND GENERATE CINDA DATA 58 

59 

1 = 1 <>0 

I READ(NIN,2018) ISTMT 61 

WRITE{LUT3) ISTMT *2 

IFtISTMT(2) .EQ. lENDX) GO TO 7 63 

1= 1*1 6A 

GO TO 1 65 

7 NSTMTS= I 66 

REWIND LUT3 67 

READINlNtlOOO) ID 68 

69 

INPUT DATA DESCRIPTION * 70 

CARD 1 FORMAT! 14, 213, 7F10.5) 71 

NSEG NO. OF SEGMENTS (MAXIMUM 30) 72 

NPAY NO. OF SURFACES, (SEGMENTS) ON PAYLOAD (MAXIMUM 10) 73 

NSPL NO. OF SHIELDS PER LAYER 74 

DELS THICKNESS OF EACH ( SHIELD+SPACER ) FT. 75 

TWALL TEMPERATURE OF TANKWALL DEG.R 76 

TSJRF TEMPERATURE OF TOP SURFACE OF SHIELD DEG.R 77 

TSPACE TEMPERATURE OF SPACE DEG.R 78 

EMISS EMISIVITY USED IN EQUATION FOR K(EFF), 79 

CONDUCTION IN NORMAL DIRECTION 80 

SPACBE EMISIVITY OF SPACE 81 

RP RADIUS OF PAYLOAD FT. 82 

83 

CARD 2 FORMAT! 6F10.5, A6) 84 

TCAROS SWITCH- TO READ INITIAL TEMPERATURES 85 

IF TCARDS =0, CARDS WILL NOT BE READ. THE PROGRAM 86 

USES TWALL AND TSURF AND DISTRIBUTES THE TEMPERATURES 87 

BETWEEN THESE TWO VALUES AS THE INITIAL TEMPERATURE 88 

DISTRIBUTION. 89 

IF TCARDS .GT. 0 THE INITIAL TEMPERATURE DISTRIBUTION 90 

IS READ FROM CAROS, (DEG. R). SEE CARDS 11 BELOW. 91 

X LENGTH OF SEMI-AXIS OF ELIPSE, ( TANKWALL ), ALONG X AXIS 92 

Y LENGTH OF SEMI-AXIS OF ELIPSE, ( TANKWALL ), ALONG Y AXIS 93 

H TANK TO PAYLOAD SPACING FT. 94 

CONI CONSTANTS USED FOR STRAIGHT LINE CURVE FIT TO COMPUTE 95 

C0N2 THERMAL CONDUCTIVITY FOR G VALUES FOR LATERAL 96 

CONDUCTION OF STANDARD NODES. K=C0N1*T+C0N2. 97 

lEXECN NAME OF EXECUTION ROUTINE CINDSS OR CINDSL 98 

CARD 3 FORMAT! 2110, F10.5) 99 

NPRINT DELTA PRINT, PRINT TEMP. DIST. EVERY NPRINT ITERATIONS 100 

NDAMP CHANG DAMPING FACTOR EVERY NDAMP ITERATIONS 101 

DELDMP DELTA TO ADD TO DAMPING FACTOR 102 

103 

CARD 4 F0RMAT(I8,9F8.3) 104 

CINCON(I) CINDA CONTROL CONSTANTS, AS REQUIRED BY EXECUTION 105 

ROUTINE BEING USED. CONSTANTS INPUT IN ORDER AS 106 

FOLLOWS ACCORDING TO FORMAT GIVEN ABOVE, 107 

NLOOP, ARLXCA, DAMPA, DRLXCA, DAMPD, TIMEO, 108 

TIMEND, OUTPUT, DTIMEI, CSGFAC. 109 

110 

CARDS 5 NTAB, ( AKVST! I ) , 1=1 ,NTAB) FORMAT! 16, 12F6. 1/ (6X, 12F6.1 ) ) 111 

TABLE OF TEMPERATURE VS. THERMAL CONDUCTIVITY FOR 112 

LATERAL CONDUCTION OF HIGH CONDUCTIVITY LAYERS. 113 

NTAB= TOTAL NO. OF ENTRIES IN TABLE (TEMPS AND K»S) 114 
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DATA FORM NTAB TEMPI Kl TEMP2 K2 TEMPN KN 115 

UNITS — K BTU/HR-FT-OEG.R 116 

T DEG.R 117 

118 

CARDS 6 ML NUMBER DF LAYERS IN EACH SEGMENT FORMAT ( NI2) 119 

120 

CARDS 7 THETAD ANGLE OF EACH SEGMENT, IN DEGREES, IF ALL SEGMENTS 121 

ARE EQUAL INPUT ONLY THE FIRST. FORMAT (N 12) 122 

THE SUM OF THE N ANGLES MUST EQUAL 90 123 

124 

CARDS 8 ON TWO CARDS (10F8.5) 125 

EMISIP(I) EMISIVITY OF EACH PAYLOAD SEGMENT (MAXIMUM 10) 126 

TPAY(I) TEMPERATURE OF EACH PAYLOAD SEGMENT (MAXIMUM 10) 127 

128 

CARDS 9 EMIST EMISSIVITY OF EACH SURFACE, RADIATION (8F10.5) 129 

130 

CARD 10 NHIK NO. OF HIGH CONDUCTIVITY LAYERS (3012) 131 

NUMHIK LAYER NUMBERS WHICH ARE HI CONO. 132 

133 

CARDS 11 NNODES IF INITIAL TEMPS READ FROM CARDS - NO. OF NODES 134 

T(I) INITIAL TEMPERATURES, DEG. R ( 1 8, 9F8. 1/10F8. 1 ) 135 

136 

READ(NIN, 1001) NSEG, NPAY,NSPL, DELS, TWALL,TSURF,TSPACE, EMI SS, 137 

1 SPACEE,RP,TCARDS,X, Y,H,C0N1,C0N2, IEXECN,NPRINT,NDAMP, 138 

2 DELDMP.CINCON 139 

DPrH=DELS*FLOAT(NSPL) 140 

THIKN= DPTH 141 

HH=H . 142 

NNSEG=NSEG 143 

NNPAY= NPAY 144 

TNODE=D.O 145 

READ(5,1003) NTAB, ( AKVST( I ) , I=1,NTAB) 146 

HIKCON=AKVST(NTAB) 147 

IF(NSEG.GT.30) NSEG=30 148 

IFMT1(2)= NBCD(NSEG) 149 

IFMT1(4)= NBCD(NSEG) 150 

READ(NIN, IFMTl) (NL( I ) , I =1 ,NSEG ) , ( THETAD ( I ) , 1=1 , NSEG ) , EMI SP,TPAY 151 

READ(NIN, 1004) { EM 1ST ( I ) , I =1 ,NSEG ) 152 

REA0(N1N, 1002) NHI K, ( NUMHIK ( I ), 1=1 , NHIK ) 153 

NSURF= NSEG + NPAY ♦ 1 154 

NSURFS=NSURF 155 

NOFA= (NSURF*(NSURF+l)/2)+2 156 

J=NSEGi-l 157 

MAXNL=0 158 

DO 201 1=1, NSEG 159 

201 IF(NL(I) .GT. MAXNL) MAXNL=NL(I) 160 

DO 200 I=J,NSURF 161 

M=I-J+1 162 

200 EMIST(I)= EMISP(M) 163 

EMIST( NSURF)=SPACEE 164 

DO 2 1*1,50 165 

2 ITEST(1)= I 166 

IF(NHIK .LE. 0) GO TO 4 167 

DO 3 1=1, NHIK 168 

M= NUMHIK(I) 169 

3 ITEST(M)= -ITEST(M) 170 

4 CONTINUE 171 

IFMTH(4)=ISPS 172 

IFdEXECN .EQ. ISL .OR. lEXECN .EQ. IWBK .OR. lEXECN .EQ. IBACK 173 

1 .OR. lEXECN .EQ. IVARB) I FMTH( 4 ) = IL PS 174 

IF(THETAD(2) .GT. 0.) GO TO 5 175 

DO 6 1=2, NSEG 176 

6 THETAD(I)= THETAO(l) 177 
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5 BTA(1)= 0. 178 

DO 9 Is2tNSEG 179 

9 BTA(I)= BTA(I-l) + THETA0(I-l» 180 

BTA(NSEG*1)»90.0 181 

WRITE(LCK2f IFMTH) 182 

WRITE(LCK2»2001) ID 183 

NSEGP1= NSEG+1 18A 

7001 NNL= I 185 

7002 NAREAL* NSEG+2 186 

7003 NAREAN* NAREAL+NSEGP I 187 

7004 NARCLM=NAREAN+NSURF+l 188 

7005 NTHETA= NARCLN+NSEGP 1 189 

7006 NITEST=NTHETA+NSEGP1 190 

7007 NKHIC= NITEST+51 191 

7008 NEMISS= NKHIC+NTAB+l 192 

7009 NAFA=NEMISS+NSURF+1 193 

7010 NSNT=NAFA+N0FA+1 194 

C 195 

C GENERATE NODE DATA 196 

C NODE NO. FOR SPACE IS NNODES 197 

C NODE NO. FOR FIRST PAYLOAD NODE IS I NNODES-NPAY ) 198 

C NODE NO. FOR TANKWALL IS (NNODES-NPAY-1 ) 199 

C 200 

C SHIELD NODES 201 

C 202 

WRITE(LCK2,2032) 203 

101 NN00ES=0 204 

205 

GENERATE INTERIOR NODES 206 

207 

DO 10 I=1.NSEG 208 

10 NNODES* NNODES + NLI I» 209 

NINSN=NNODES 210 

211 

ADD NO. OF SURFACE NODES TO NNODESt SURF NODES ONLY ALONG TANKWALL 212 

213 

NNODES»NNODES+NSEG 214 

TEMP* TNODE 215 

WRITE! LCK2,2003»M1 .NNODES. Ml . TEMP. ONE . ONE. ONE . ONEM 216 

217 

TANKWALL NODE 218 

219 

NNODES* -NNODES-1 220 

TEMP* TWALL 221 

WR1TE(LCK2.2004) NNODES. TEMP .ONE 222 

223 

PAYLOAD NODES (BOUNDARY NODES) 224 

225 

DO 105 1*1.NPAY 226 

NN0DES*NN0DES-1 227 

TEMP*TPAY(I) 228 

WRITE(LCK2.2004) NNODES. TEMP .ONE 229 

105 CONTINUE 230 

231 

NODE FOR SPACE (BOUNDARY NODE) 232 

233 

NNODES = NNODES-1 234 

TEMP* TSPACE 235 

WRITE(LCK2.2004) NNODES.TEMP.ONE 236 

104 WRITE(LCK2.2007) 237 

NNODES* lABS(NNODES) 238 

IF(TCARDS .EQ. 0.) GO TO 115 239 
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READC MIN, 1005) NNSPCH, (T ( I) , 1=1, NNSPCH ) 240 

GO TO 110 241 

115 CONTINUE 242 

N=0 243 

IF(TSJRF) 207,207,208 244 

245 

CONSTANT INITIAL TEMPS 246 

247 

207 TS=ABS(TSURF) 248 

DELT=0.0 249 

GO TO 205 250 

251 

DISTRIBUTE TEMPS BETWEEN TWALL AND TSURF 252 

253 

208 DEL T=( TSURF-TWALL ) /FLOAT (MAXNL ) 254 

TS=TWALL 255 

205 D0210 I=1,NSEG 256 

NLAY=NL(IH-1 257 

N=N^-1 258 

T(N)=TS 259 

D021I J=2,NLAY 260 

N=N+1 261 

T(N)=TCN-1) +DELT 262 

211 CONTINUE 263 

210 CONTINUE 264 

DO 212 1=1, N 265 

212 T( I)=T( I )-TABS 266 

110 T(NNODES)=TSPACE -TABS 267 

M=NN0DES-NPAY-1 268 

DO 213 I=1,NPAY 269 

J=M^*I 270 

213 T( J)=TPAYC n-TABS 271 

T(M)=TWALL-TABS 272 

WRITECLUT5) ( T ( I ) , 1= 1 , NNODES ) 273 

C 274 

C COMPUTE AREAS FOR LATERAL CHORIZ), AND NORMAL (VERT), CONDUCTORS 275 

C PUT AREAS INTO ARRAYS AREAL, AND AREAN AND OUTPUT IN ARRAY DATA AS 276 

C ARRAY 2 AND ARRAY 3 RESPECTIVELY. 277 

C ARRAY ARCLN, LENGTHS OF NODES, IS ARRAY 4. 278 

C 279 

80 CONTINUE 280 

ASQ=X*X 281 

BSQ=Y*Y 282 

ASQBSQ= ASQ*BSQ 283 

ERAD(1)=Y 284 

K=0 285 

DO 20 I=I,NSEG 286 

ANG1= BTA( I) 287 

ANG2= 8TA( I+l) 288 

THETAt I)= THETADI I )*0. 01745329 289 

M= THETAD(I) 290 

DLAY=NL( I ) 291 

DO 81 L=1,M 292 

K=K+1 293 

81 SLAY(K)=DLAY 294 

ANG2= ANG2*0. 01745329 295 

SINAL= SIN(ANG2) 296 

COSAL= C0S(ANG2) 297 

SSO= SINAL**2 298 

CSQ= CDSAL**2 299 

ERAD(M^1)= SQRT(ASQBSO / (BSO*SSO ♦ ASQ*CS0) ) 300 

DMLI= DPTH*FL0AT(NL( I ) ) 301 

Rl= ERAD( I )+DMLI 302 

R2= ERAD( I+1)+0MLI 303 
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RAD= (Rl+R2)/2* 304 

ARCLN( I )= RAD*THETA( I ) 305 

SMALLR= R2* SINAL 306 

ENDW= SMALLR* PI04 307 

AR£AL(n= ENDW • DPTH 308 

20 CONTIf^UE 309 

CALL VIEWF 310 

311 

GENERATE CONDUCTOR DATA , ALL CONDUCTANCE VALUES WILL BE COMPUTED 312 

IN VARIABLES 1 (THEREFORE 1.0 WILL BE USED 313 

IN THIS ROUTINE FOR C. ) 314 

315 

WRITE(LCK2,2005) 316 

NORMAL CONDUCTORS IN SHIELD 317 

318 

N0DA=1 319 

NC0ND=1 320 

DO 11 I=ltNSEG 321 

NLAY=NL( I ) 322 

N0DB=N0DA+1 323 

WRITE(LCK2»2006)NC0ND,NLAY,M1,N0DA,M1»N0DB.M1,0NE,0NE*0NE.0NE 324 

N0DA=N0DA+NLAY+1 325 

11 NCOND=NICOND+NLAY 326 

327 

LATERAL CONDUCTORS IN SHIELD 328 

329 

N0DA=2 330 

DO 15 I=2,NSEG 331 

NLAY=NL(I) 332 

NLM1P2=NL ( I-l)^l 333 

N0DB=N3DA+NLM1P2 334 

WRITE( LCK2,2006)NCOND,NLAY,M1,NDDA,M1,NODB,M1,ONE,ONE,ONE,ONE 335 

NODA=NODB 336 

15 NCOND=NCOND+NLAY 337 

338 

CONDUCTORS FROM EDGE OF SHIELD TO TANKWALL 339 

340 

NODA=0 341 

N0DB=NN0DES-NPAY-1 342 

NCOND=NCOND-1 343 

144 DO 14 I=ltNSEG 344 

N0DA=N0DA+1 345 

NC0ND=NC0ND+1 346 

WRITE( LCK2,2026) NCDND , NODA, NODB, TEN20 347 

14 NODA= NOOA+NLd) 348 

349 

COMPUTE TABLE OF SURFACE NODE NUMBERS 350 

351 

NODA=D 352 

DO 16 I=1,NSEG 353 

NODA= NODA+NLI I )+l 354 

NSNODd )= NODA 355 

16 CONTINUE 356 

J=NSEG+1 357 

JJ=NNODES*NPAY~J 358 

DO 162 I=JtNSURF 359 

162 NSNODt I )=JJ+I 360 

WRITE(LCK2.2007) 361 

362 

GENERATE CONSTANTS DATA 363 

364 

WRITE! LCK2,2008) ( I « C d ) « I ^ 1 » 20 ) 365 
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WRITE(LCK2,2024) ( I 9 C ( I ) t I=2I» 28 ) 366 

WRITE(LCK2f2010) ( I» C ( I ) 9 1=29, 36) 367 

WRITE(LCK2,2022) ( CONNAM ( I ) , CINCONC Dtl^IflO) 368 

WRITECLCK2, 20071 369 

C 370 

C GENERATE ARRAY DATA BLOCK 371 

C TO OUTPUT ARRAY CALL SUBROUTINE RI TAR Y (N • A ,N0# NF , ZN) 372 

C WHERE N=ARRAY NUMBER 373 

C A=ARRAY TO BE PUT OUT 374 

C N0= NUMBER OF ELEMENTS IN THE ARRAY 375 

C NF= EITHER 1 OR 2 NF*1 IF ARRAY IS FLOATI 376 

C NF=2 IF ARRAY IS INTEGE 377 

C ZN=BCD ARRAY WHICH BECOMES A TITLE FOR 378 

C ARRAY N. 379 

C 380 

WRITE(LCK2,2009) 381 

CALL RITARY(lfNL,NSEG,2,21) 382 

CALL RITARYC2, AREAL, NSEG»I,Z2) 363 

CALL RITARY(3,A,NSURF,1,23) 384 

CALL RITARY(4,ARCLNf NSEG,1,Z4) 385 

CALL RITARYC 5,THETAD,NSEG, 1,Z5) 386 

CALL RITARY( 6 ,ITEST, 50 , 2 ,Z 6 ) 387 

CALL RITARY(7,AKVSTf NTAB,1,Z7) 388 

CALL RITARY( 8 ,EMIST,NSURF, l,Z 8 ) 389 

CALL RITARY ( 9 , AFA,NOFA, 1, Z9 ) 390 

CALL RITARY(10,NSNOD,NSURF,2,Z10) 391 

WRITE(LCK2,2007) 392 

C 393 

C GENERATE EXECUTION BLOCK 394 

C 395 

WRITE (LCK2,20I2) 396 

WRITE(LCK2,2020) 397 

WRITE(LCK2,2033) lEXECN 398 

WRITECLCK2,2007) 399 

C 400 

C GENERATE VARIABLES 1 BLOCK 401 

C 402 

WRITE(LCK2,2015) 403 

WRITE(LCK2,2007) 404 

C 405 

C GENERATE VARIABLES 2 BLOCK 406 

C 407 

WRITE(LCK2,2014) 408 

WRITE(LCK2,2031) 409 

WRITE(LCK2,2007) 410 

C 411 

C GENERATE OUTPUT CALLS 412 

C 413 

WRITE(LCK2,2016) 414 

WRITE(LCK2,2029) 415 

WRITE(LCK2,2028) 416 

WRITE(LCK2,2025) 417 

WRITE(LCK2,2030) 418 

WRITE(LCK2,2007) 419 

C 420 

C WRITE SUPPLIED SUBROUTINES ONTO TAPE LCK2 421 

C 422 

WRITE! LCK2,2017) 423 

DO 60 I=1,NSTMTS 424 

READ(LUT3) ISTMT 425 

60 WRITE(LCK2,2018) ISTMT 426 

WRITE(LCK2,2019) 427 

REWIND 2 428 
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429 

WRITE SWITCH SYSINl TO SYSUT6, TO READ CINDA DATA BY CINDA 430 

431 

END FILE LCK2 432 

WRITE(LCK2,2021) 433 

END FILE LCK2 434 

END FILE LCK2 435 

END FILE LCK2 436 

REWIND LCK2 437 

WRITE! NOUTtSOOO) 438 

5000 F0RMAT(51H END OF PREPRE PGM. NOW SWITCHTO INPUT FOR PREPRO) 439 

WRITE! 6t5001) NSTMTS 440 

5001 F0RMAT!45H NO. OF STATEMENTS IN SUBROUTINES LOADED = 15) 441 

CALL EKIT 442 

READ!NINtlOOO) ID 443 

1000 FORMAT! 13A6.A2) 444 

1001 FORMAT!I4«2I3»7F10.5/6F10.5,A6/2I10.F10.5/I8«9F8.3) 445 

1002 F0RMAT!30I2) 446 

1003 F0RMAT!I6»12F6.1/!6X»12F6.1) I 447 

1004 FORMAT! 8F10. 5) 448 

1005 FORMAT! 18, 9F8.1/!10F8.1) ) 449 

2000 FORMAT! 20HSEXECUTE IBJOB, 64X / 450 

1 34H$ID C0W6ILL PREPRO, 50X/ 6H$IBJ0B ,78X / 451 

2 27H$IEDIT SYSCK1,SCHF2 ,57X / 13HSIBLDR PREPRO ,71X / 452 

2 13H$IBLDR TAPDEF,71X/ 13H$IBLDR FVIO. ,71X/ 453 

3 27H$DRIGIN ALPHA, SYSLB2 ,57X / 13H$IBLDR GENLNK ,71X / 454 

4 27H$DRIGIN ALPHA, SYSLB2 ,57X / 13HSIBLDR PSEUDO ,71X / 455 

5 13H$IBL0R PACK43 , 71X / 13HSIBLDR ORMIN , 71X / 456 

6 27HJDRIGIN ALPHA, SYSLB2 ,57X / 457 

7 13H$IBLDR CODERD , 71X / 13H$IBLDR OATARD , 71X / 458 

8 27H*DIUGIN ALPHA, SYSLB2 ,57X / 13H*IBLDR PRESUB ,71X / 459 

9 13HSIBLDR CIN0A4 , 71X / 13H$IBLDR SEARCH, 71X / 460 

A 27H$0RIGIN ALPHA, SYSLB2 ,57X / 13HSIBLDR COPYSB, 71X / 461 

B I3H$IBLDR LOCATE , 71X / 13H$IBL0R FIBLOR ‘, 71X / 462 

C 13H$IBLDR READTP , 71X / 13HSIBLDR CIN90 , 71X / 463 

D 27HSDRIGIN ALPHA, SYSLB2 ,57X / 13HSIBLDR INITAL , 71X / 464 

E 27H$DRIGIN ALPHA, SYSLB2 ,57X / 13HSIBL0R FINAL , 71X / 465 

F 27HS0RIGIN ALPHA, SYSLB2 ,57X / 466 

G 13H$IBLDR SKIP , 71X / 13HtIBLDR SPLIT , 71X / 467 

H 5HSDATA , 79X / 80X ) 468 

2001 F0RMAT!7X,5HBCD 9,9A6, 14X/7X,5HBCD 5, 5A6 , 38X/7X , 3HEND, 70X ) 469 

2003 FORMAT! 470 

1 7X,4HGEN ,I1,1H,,I4,1H,,I1,1H,,F7.2,1H,,F3.1,1H,,F3.1,1H,,F3.1, 471 

2 IH, ,F4.1,36X) 472 

2004 F0RMAT!11X,I5,1H,,F6.1,1H,,F4.1,52X) 473 

2005 FORMAT! 7X,19HBCD 3CDNDUCT0R DATA ,54X) 474 

2006 F0RMAT!7X,4HGEN , 7 ! I 4, IH , ) , 3 ! F3. 1 , IH, ) , F3. 1 , 19X) 475 

2007 F0RMAT!7X,3HEND,70X) 476 

2008 F0RMAT!7X,19HBCD 3C0NSTANTS DATA,54X / 477 

1 !11X,4!I3,1H,, 15,1H, , ), I3,1H, , I5,20X) ) 478 

2009 F0RMAT!7X,15HBCD 3ARRAY DATA,58X) 479 

2010 FORMAT! !11X,3! I2,1H, ,F13.5,1H,, ),I2,1H,,F13.5,2X) ) 480 

2012 FORMAT ! 7X, 14HBCD 3EX ECUT ION, 59X/23HF DIMENSION X!3000),57X / 481 

1 15HF NDIM=3000,65X / 12HF NTH=0,68X ) 482 

2014 FORMAT !7X,16HBCD 3VAR1ABLES 2,57X) 483 

2015 F0RMAT!7X,16HBCD 3VARIABLES 1, 57X /15HF CALL GETG,65X) 484 

2016 F0RMAT!7X,17HBCD 30UTPUT CALLS, 56X) 485 

2017 F0RMAT!7X,16HBCD 3SJBR0UT1NES ,57X ) 486 

2018 FORMAT! 13A6,A2) 487 

2019 FORMAT ! 7X, 16HBC0 3EMD OF DATA ,57X) 488 

2020 F0RMAT!11HF NN=10,69X/22HF IF ! NN.GT. 0) 60T03, 58X / 489 

1 11X,24HLAGRAN!RTEST,STEST,A7,7) ,45X/ 490 
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2 llX,22HrRRADI(A3,A8,A9,STEST» ,47X / 491 

2 14HF 3 CONTINUE, 66X / 492 

2 16HF CALL TDIST.64X) 493 

2021 F0RMAT(6H$IBSYS,74X/21H$REWIN0 SYSLB3.59X/ 494 

1 21H$REWIND SYSUT6, 59X/28H$SWITCH SYSIN1.SYSLB3.52X/ 495 

2 28H$SWITCH SYSIN1,SYSCK2,52X) 496 

2022 F0RMAT(11X,A6,1H,,I4,3(1H,,A6,1H,,F9.6),7X / 497 

1 11X,A5,1H,,F9.6,3(1H,,A6, 1H,,F9.4),2X / 498 

2 11X,A5, IH, ,F9. 6, IH, ,A6, 1H..F9. 6,36X1 499 

2023 FORMAT! 12X,3HEND,65X) 500 

2024 FORMAT! (11X,3!I2,1H,,1PE13.6,1H,,),I2,1H,,1PE13.6,2X)I 501 

2025 F0RMAT!26HF WRITE !6, 1000) L00PCT.54X / 502 

1 28HF1000 F0RMATI9HL LOOPCT= I6),52X) 503 

2026 F0RMAT!11X,3!I5,1H,I,3X,A6,42X) 504 

2027 F0RMAT!11X,3!I5,1H,),E16.8,35X) 505 

2028 FORMAT! 11X,6HPRNTMP,63X) 506 

2029 FORMAT! 13HF NN=K!4),67X/ 507 

1 18HF DO 10 1=1, NN, 62X / 508 

2 22HF 10 T!I)=T!I )*460.0 ,58X) 509 

2030 FORMAT! 18HF DO 11 I=1,NN,62X/ 510 

2 22HF 11 T! I)=T!I )-460.0 ,58X) 511 

2031 FORMAT (17HF CALL EN00UT.63X) 512 

2032 FORMAT! 7X, 14HBCD 3N30E DATA,59X) 513 

2033 FORMAT! 11X,A6,63X) 514 

END 515 


$IBFTC VIEWFC 


SUBROUTINE VIEWF 1 

RADIATION YIEW FACTOR BETWEEN FLAT SURFACED PAYLOAD AND SPHEROIDAL 2 

TANK WITH VARYING ML I THICKNESS. 3 

INPUT REQUIRED 4 

SLAY - NO. OF NODES IN EACH 1 DEGREE SEGMENT 5 

RP - PAYLOAD RADIUS 6 

X AND Y - MAJOR AND MINOR RADII OF TANK IWITHOUT COVERING) 7 

NSEG - NO. OF SEGMENTS ON TANK 8 

NPAY - NO. OF RINGS ON PAYLOAD 9 

NSURF- TOTAL NO. OF SURFACES IN RADIATION ENCLOSURE 10 

H - TANK TO PAYLOAD SPACING 11 

COMMON /RADIO/ NSEG, NPAY, NSURF , X,Y,H,RP, THIKN, PI , SLAY, THETA, A, AF A 12 

COMMON /TAPE/ NIN,N0UT,LCK2 13 

DIMENSION SLAY!90),A!41) ,AFA!1681),THETA!30) 14 

DIMENSION SAREA!90),DAREA!90),SEGA!100),DF!600),FI600),R!10) 15 

DIMENSION SAF190,100),AFI41,41),RVF!41,41) 16 

DATA H1,W2,W3,W4,W5,W6,EY1,EY2 /6H A,6H R,6H DAREA, 17 

1 6H SEGA,6H RFV,6H AF,6H!I) ,6H!I,J) / 18 

T0PI=2.*PI 19 

ONED= Pl/180. 20 

HALFD=0NED/2. 21 

RPSQ=RP*RP 22 

N=10.*RP 23 

PSI=-HALFD 24 

DO 28 K=l,90 25 

PSI=PSI*ONEO 26 

10 DMLI= SLAY(K)*THIKN 27 

AX=X+DMLI 28 

BX=Y+DMLI 29 

AXSQ=AX»AX 30 
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BXSQ=BX*BX 
SPSI = SIN(PSI ) 

CPSI=CDSCPSI ) 

SSQ=SPSI*SPSI 

CSQ=CPSI*CPSI 

AXBX=AX*BX 

Rl= AXBX*SQRT(1./(BXSQ*SSQ+AXSQ*CSQ) ) 

Y1«R1*SPSI 

Z1=R1*CPSI 

DH=BX-Z1 

XM1=(AXSQ*ZU/(BXSQ*Y1) 

WSQ=(AXSQ-BXSQ)/BXSO/BXSQ 

W=SQRrCWSQ) 

H0W=0.5/W 

SRQB=SQRTC 1,+WSQ*BXSQ) 

RJ=BX/2.*SRQB 
S= HDH*ALOG(BX*W'^SRQB) 

UPLIM=RJ+S 

PSIX^PSH-HALFD 

SPX=SIN(PSIX) 

CPX^COSCPSIX) 

SPXSQ=SPX*SPX 

CPXSQ=CPX*CPX 

R1X= AXBX*SQRT(1./(BXSQ*SPXSQ+AXSQ*CPXSQ) ) 
Z1X=R1X*CPX 

SRQZX=SQRT(WSQ*Z1X**2+1. ) 

T=Z1X/2.*SRQZX 
U= HOW»ALOG(ZlX*W+SRQZX) 

BOTLIM= T+U 

SAREA( X)=TOPI*AX*(UPL IM-BOTL IM) 

IF(K.EQ.l) GO TO 12 

DAREA( K)=SAREA(K )-SAREA(K-U 

GO TO 14 

12 OAREAU) = SAREA(K) 

14 DO 28 1 = 1, N 
RING=I 

RHO=.l*(RING-*5) 

SIXRNG=6.*RING 

SEGA (I )=PI*( t .l*RING)**2-( .1*(RING-1.) )**2)/SIXRNG 

NSGT=S*I 

GAMA=0. 

DO 24 J=1,NSGT 

IF( J.EQ.l ) GO TO 16 

GAMA=GAMA+(TOPl/SIXRNG) 

16 X2=RHD»COS(GAMA) 

Y2=RH3*SIN<GAMA) 

XA= X2 

YB= Y2-Y1 

ZC= H-DMLI+DH 

RV=SQRT(XA**2+YB**2+ZC**2) 

SM= SQRTt 1,+XM1**2) 

COSPl= ( YB+XMl*ZC)/( RV*SM) 

C0SP2= ZC/RV 

IF(COSPl*LE.O. ) GO TO 18 

OF( J)=(COSPl*COSP2*SEGAC I) )/(PI*RV*RV) 

GO TO 20 
18 DF(J)=3. 

20 IF{J.EQ*l) GO TO 22 
F(J)=0FCJ)+F(J-1) 

GO TO 24 
22 F(J)=DF(J) 

24 CONTINUE 

SAF(K, I ) = DAREACK)*FI NSGT ) 
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28 CONTIMJE 94 

XNPAY-MPAY 95 

DO 30 II=1,NPAY 96 

XI=II 97 

RR=SQRT( i RPSQ*XI l/XMPAY) 98 

NDELR=10.*RR 99 

XNOELR=NDELR 100 

R(in = .l*XNOELR 101 

JJ=II + ^ISEG 102 

IF(II.EQ.l) GO TO 29 103 

A( JJ)= PI*(R(II)**2-R( II-1)**2) 104 

GO TO 30 105 

29 A( JJ)=PI*RCII)**2 106 

30 CONTIJ^JE 107 

DO 304 KK=ltNSEG 108 

IF(KK.EQ.l) GO TO 301 109 

NDEG1=MDEG2+1 110 

NDEG=THETA(KK) 111 

NDEG2==MDEG1+NDEG-1 112 

GO TO 302 113 

301 N0EG1=1 114 

NDEG2=THETA(KK) 115 

302 ASUM=0, 116 

DO 303 K=NDEG1,NDEG2 117 

ASUM=ASUM+DAREA(K) 118 

303 C0NTI*^IJE 119 

A(KK)=ASUM 120 

304 CONTIMJE 121 

A(NSURF)= 122.211 122 

DO 40 XK=1,NSEG 123 

IF(KK.EQ.n GO TO 32 124 

NDEG1=MDEG2+1 125 

NDEG=THETA(KK) 126 

NDEG2=MDEG1+NDEG-1 127 

GO TO 34 128 

32 NDEG1=1 129 

NDEG2=rHETA(KK) 130 

34 DO 40 II=1,NPAY 131 

IF(II.EQ.l) GO TO 36 132 

NDR1=MDR2+1 133 

NDR2=1D.*RCII) 134 

GO TO 37 135 

36 NDR1=1 136 

NDR2 = 1D.*R( I I ) 137 

37 SUM=0. 138 

DO 38 K=NDEG1»NDEG2 139 

DO 38 I=NDR1*NDR2 140 

SUM=SJM+SAF(K,I) 141 

38 CONTIMJE 142 

JJ=II4^MSEG 143 

AF(KK, JJ)=SUM 144 

40 CONTIMJE 145 

DO 42 I=1,NSEG 146 

DO 42 J=ItNSEG 147 

42 AF(I,J)=0. 148 

NS1=:NSEG+1 149 

NS2=NSJRF-1 150 

DO 44 I=NS1.NS2 151 

DO 44 J=I,NS2 152 

44 AFCItJ)=0. 153 

AF(NSJRF,NSURF)=0. 154 

DO 48 I^ltNSEG 155 

SUM=0 156 
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46 

48 


50 

52 


60 


100 


200 


1100 

1001 

1002 


DO 46 K=NS1»NS2 
SUM=SJ'<+AF(I,K) 
AFCI*^SURF)=A(I)-SUM 
DO 52 I=NS1,NS2 
SUM=0 

DO 50 <=lfNSEG 
SUM=SJM+AF(K.I) 

AF(I.'^SURF) = A{I)-SUM 
DO 60 KK=1,NSEG 
DO 60 JJ=NSltNSURF 
RVF(K<, JJ)=AFIKK, JJ)/A(KK) 
AFA(1)=0.0 
AFA(2)=0.0 
1 = 2 


DO 100 KK=1,NSURF 
DO 100 JJ=KK,NSURF 
1= I+l 

AFAU)= AFtKKfJJ) 
Z1 = W1 


Z2=EY1 

WRITE(MOUT,1001) 
Z1 = W2 

WRITE(MOUT,1001) 
Z1 = W3 

WRITEt ^OUT.lOOl) 
Z1 = W4 

WRITE(M0UT,1001) 

Z1 = W5 

Z2=EY2 

WRITE( ^J0UT,1001) 
DO 200 KK=1,NSEG 
WRITE(M0UT,1002) 
CONTIMJE 


Zl«Z2t(A( n ,I=1,NSURF) 
Z1.Z2»(RI D*I = ltNPAY) 
Z1,Z2, C DAREA( I) , 1 = 1,90) 
Zl,Z2f ( SEGA( I )t I=1,N) 

Zl, Z2 

(RVF(KK,JJ) ,JJ=NS1,NSURF) 


Z1 = W6 

WRITE(>40UT,1001) Z1,Z2 
DOllOO KK=1,NSURF 

WR I TE( MOOT , 1002 ) ( AF ( KK, JJ ) , J J=KK, NSURF ) 

CONTIMUE 

RETURN! 

FORMAT ( lHK,2A6/( IX, 1P10E13.5) ) 

FORMAT! IX, 1P10E13.5) 

END 
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$IBFTC RAYOJT 


SUBROUTINE R IT AR Y ( NA , A , N , J , HE AD ) 1 

DIMENSION Ad), IFMTK5), IFMT2(7) 2 

DIMENSION HEA0(9) 3 

COMMON/BCD/ NBCD(60) 4 

COMMON /TAPE/ N IN, NOUT , LCK2 5 

DATA IFMTl(l) /30H(12X,4 tlH,), 4X ) / 6 

DATA IFMT2d) /42H(12X, X) / 7 

DATA IFLOTl , IFLDT2 , IFL0T3 8 

1 / 6H(F12.7 , 6HF12.7, , 6H, IH. ) , / 9 

DATA IFIXl , IFIX2 , IBLANK 10 

1 / 6H (112 t 6H 112, , 6H / 11 

DATA lEl • IE2 12 
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6HE12.5 


1 / 6H(E12.5 • 6HE12.5, / 

NSEG= M 
JJ= J 

IFMT2(3)= IBLANK 
IFMT2(<^)= IBLANK 
ITEM2= IFL0T3 
GO TO (10, 11.12), JJ 
12 IFMT1(2)=IE1 
IFMTl(^f)=IE2 
IFMT2(5)=IE2 
ITEM1=IE1 
GO TO 21 

10 IFMT1(2)= IFLOTl 
IFMT1(4)= IFL0T2 
IFMT2(5)= IFL0T2 
ITEM1= IFLOTl 

GO TO 21 

11 IFMT1(2)= IFIXl 
IFMTl(^fr)= IFIX2 
IFMT2(5)= IFIX2 
ITEM1= IFIXl 

21 Nl= 1 
N2= 5 

WRITEi LCK2,2011) NA, ( HEAD( I) , I =* 1 , 9 ) 
31 IF (N2 .GE. NSEG) GO TO 40 

WRITE (LCK2,IFMT1) ( A ( I) , I =N 1 , N2 ) 
Nl= N1^5 
N2= N2^5 
GO TO 31 
40 NN= NSEG-Nl 
N2= NSEG 

IF(NN .LE. 0) GO TO 39 
IFMT2(3)* ITEMl 
IFMT2(4)= ITEM2 
IFMT2(2)= NBCD(NN) 

NBL= 56-13*NN 
42 IFMT2(6)= NBCDCNBL) 

WRITE(LCK2,IFMT2) ( A ( I ) , I = N1,N2 ) 
WRITE (LCK2,2013> 

RETURN 


39 NN=1 

IFMT2(2)=NBCD(NN) 

NBL=56 
GO TO 42 

2011 FORMATt IIX, I2.3H S,9A6,10X) 

2013 F0RMAT(12X,3HEN0,65X) 

END 
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SIBFTC COMDAT 


1 

2 

3 

4 

5 


BLOCK DATA 


COMMO'J 

/BCD/ 

NBCD(60) 


DATA NBCD( 1) 

/360H 

1 
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12 
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END 
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58 

59 

60 


1 

2 


6 

7 

8 

3 

17 
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SIBFTC CNDCTR 


SUBROUTINE GETG 1 

COMMON /KONST/ NSEGtNPAY.NSURF,NNOOEStMAXNL.NPRINTtNOAMP,NNL» 2 

1 NAREAL*NAREAN*NARCLN.NTHETA,NITEST.NKHICtNEMISStNAFA(NSNTf 3 

2 NEXTRA(3)t ^ 

3 S IGMA. OPTH. DELS » EMI SStCQNl* CON2tOELOMP,X« TABS t HI KCONtNPSWCH, RAD, H 5 

COMMON /CONO/ G(l) 6 

COMMON /TEMP/ TU) 7 

COMMON /SOURCE/ Q(l) 8 

COMMON /ARRAY/ All) 9 

COMMON /FIXCON/ KAfl) 10 

COMMON/QCOM/ QTS(30) , QN0RM(30) 11 

DIMENSION SAREA(41).EMISIV(41).T0(41),NT0(1).NARYI1I 12 

EQUIVALENCE (A, NARY), (TQ,NTQ) 13 

EQUIVALENCE (NLOOP,KA( 5 ) ) , ( OAHPA.KAC 9) ) , (LOOPCT, KA( 20) ) 14 

EQUIVALENCE ( ARLXCC, KA ( 30) ) 15 

DATA NTIME /!/ 16 

NPSWCH=0 17 

IF(LOOPCT.GE.NLOOP) CALLPCHT(l) 18 

IFINPSWCH .GT. 0) GO TO 9998 19 

IFIARLXCC .LT. .01) DAMP A=DAMPA+DELDMP 20 

IFIDAMPA .GT. 1.0) DAMPA=1.0 21 

IFIMOD(LOOPCT,10O ) .NE. 0 ) GO TO 32 22 

CALL TIMLFT(TLEFT) 23 

TLEFTM=TLEFT/3600. 24 

WR1TEI6,9999) TLEFTM, LOOPCT 25 

IFITLEFT .LT. 14400.0) CALL PCHTI2) 26 

IFINPSWCH .GT. 0) GO TO 9998 27 

32 CONTINUE 28 

IFIMOD(LOOPCT,NPRINT) .NE. 0) GO TO 33 29 

35 NSP1=NSURF+1 30 

CALL ENDOUT 31 

WRITE<6,1001) LOOPCT, DAMPA,ITQtI), 1=1, NSPl) 32 

33 CONTINUE 33 

NRAD=RAD 34 

35 

GET G FOR ALL NORMAL (VERTICAL) CONNECTORS 36 

37 

AKPART= OELS*SIGMA/( (2./EMISS)-!.) 38 

NODA= 3 39 

NCONO= 0 40 

Nl= 1 41 

DO 11 I=1,NSEG 42 

M=NNL*I 43 

NLAY= NARY(M) 44 

MA= NAREAN+I 45 

AOD=A( MA) /DPTH /8. 46 

N2 = NH-NLAY-1 47 

30 DO 12 J=N1,N2 48 

NOOA= NODA+1 49 

NODB= NODA+l 50 

Tl= T(NODA) +TABS 51 

T2= T(NODB) -^TABS 52 

TPART= (T1+T2)*(T1*TI*T2*T2) 53 

AK= AKPART*TPART 54 

G(J)=ADD*AK 55 

12 CONTINUE 56 

NT2=N0DA 57 

41 N3=N0DB-NLAY 58 

N4=N3+1 59 

G(N1)=G(N1)*2. 60 

QNORHl I )=G(N1)*(T(N4)-T(N3) ) 61 
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43 CONTINUE 62 

Nl= 63 

40 TSURF=T(NT2) 64 

11 NODA= NODA+1 65 

C 66 

C COMPUTE G FOR LATERAL CONNECTORS 67 

C 68 

NODA= 1 69 

NCONO= N2 70 

DO 15 I=2,NSEG 71 

M=NNL4-I 72 

NLAY= NARY(M) 73 

NLM1P2=NARY(M^1)+1 74 

L=NARCLN+I 75 

ELl=A{L-l)/2. 76 

EL2=A(L)/2. 77 

LL=NAREAL+I-1 78 

AREA=A(LL) 79 

DO 13 J=1,NLAY 80 

NODA= NODA+l 81 

NODB= NODA+NLM1P2 82 

NCOND= NCOND^-l 83 

Tl= T(NODA) +TABS 84 

T2= T(NODB) +TABS 85 

PM=NITEST+J 86 

ITEST=NARY(MM) 87 

IF(ITEST)20, 19,19 80 

20 PMM=NKHIC 89 

IF(T1 .GT. 240.) GO TO 18 90 

CALL LAGRANCTi,AKl,AIMMM),7) 91 

17 IF(T2 .GT. 240.) GO TO 16 92 

CALL LAGRAN{T2,AK2,A(MMM),7) 93 

GO TO 21 94 

18 AK1=HIKC0N 95 

GO TO 17 96 

16 AK2 = HIICC0N 97 

GO TO 21 98 

19 AK1 = CDNI1*T1 + C0N2 99 

AK2=C0N1*T2 + C0N2 100 

21 DEN= ELl/AKl + EL2/AK2 101 

G(NCOND)= AREA/DEN 102 

13 CONTINUE 103 

15 N0DA=NDDA+1 104 

105 

COMPUTE SOURCE TERMS FOR NET Q, RADIOSITY, USING SUBR. IRRAOI FOR 106 

INTERCHANSE BETWEEN FACE OF SHIELD AND PAYLOAD 107 

108 

GO TO (600,601),NTIME 109 

600 NTIME=2 110 

J=NSURF+1 111 

M=NAREAN-1 112 

K=NEMISS-1 113 

DO 602 1 = 1, J 114 

PM=M+I 115 

KK=K+I 116 

SAREA( I )=A(MM) 117 

EMISI^( I )=A(KK) 118 

602 CONTINUE 119 

NARY(NAFA+1)=NSURF 120 

A(NAFA^-2)=SIGMA 121 

WRITE(5,1000) (SAREAI I ), 1=1, J) 122 

WRITE(6,1000) (EMISIV(I),I=1,J) 123 
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601 CONTINUE 124 

DO 500 I=1,NSURF 125 

K=NSNTH 126 

J=NARY(K) 127 

500 TQn+ll=T(d) 128 

NTQ(1)=NSURF 129 

CALL IRRAOI(SAREA»EMISIV*A(NAFA)tTQI 130 

00 505 I=1,NSEG 131 

K=NSNT+I 132 

J=NARY(K) 133 

0(J)=Q( J)+TQ(l-i^l) /8. 134 

QTS(I)=Q(J) 135 

505 CONTINUE 136 

9998 RETURN 137 

1000 FORMAT! IX, 10G13.5) 138 

1001 FORMAT! lOHL LOOPCT= 16, 5X,6HDAMPA= F6.4, 139 

1 /22H NET Q FROM RAD lOS ITY/ ! IX, 10G13. 5 ) ) 140 

9999 FORMAT! 12H TIME LEFT= F15.2, lOH LOOPCT= 18) 141 

END 142 


$IBFTC OUTENO 

SUBROUTINE ENDOUT 1 

COMMON/TITLE/ HEAOI20) 2 

COMMON /KONST/ KID 3 

COMMON/FIXCON/ X!27).L,NP 4 

COMMON/QCOM/ QTS 1 30) , QNORM ! 30 ) 5 

DIMENSION 0TS8!30),QN8!30) 6 

NSEG=K!1) 7 

STS=0. 8 

SNO=0. 9 

CALL PICTR 10 

DO 10 I=1,NSEG 11 

0TS8!I )= QTS! I )*8. 12 

QN8! I) = QNORM! I )«8. 13 

STS=STS+QTS! I ) 14 

10 SNO=SND+QNORM! I ) 15 

STS8=STS*8. 16 

SN08=SN0*8. 17 

WRITE!6,2000) !HEA0! I ) , I =1, 20) , ! I , QTS! I ) , QNORM ! I ) , 18 

1 QTS8!I ),QN8!I),I=1,NSEG) 19 

WRITE1S,2001) STS,SN0,STS8,SN08 20 

N=K!31) 21 

L=0 22 

NP=0 23 

IF!N ,E0. 0) GO TO 4 24 

GO TO !1.2),N 25 

1 WRITEI6,2003) 26 

GO TO 3 27 

2 WRITEI6,2002) 28 

3 WRITE!6,2004) 29 

4 RETURN 30 

2000 FORMAT! 1H1.10X,20A6/9HL NO. SEG,8X,3HQTS,11X,5HQN0RM, 31 

1 24X,5HQTS«8,9X,7HQN0RM«8 // 32 

2 I 18,4X.2F14.8,15X,2F14.8) ) 33 

2001 FORMAT! 12HLSUMMATI0N , 2G14. 5, 15X,2G14. 5) 34 

2002 FORMAT! 71HLN0TE . . . . THIS CASE STOPPED FOR TIME EXCEEDED, IT HAS NO 35 

IT YET CONVERGED. ) 36 

2003 FORMAT!52HL....THIS CASE STOPPED FOR ITERATION COUNT EXCEEDED. ) 37 

2004 FORMATI53HK ....CAROS WERE PUNCHED TO RESTART FROM THIS POINT. ) 38 

END 39 
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SIBFTC DISTRT 


SUBROJTINE TOIST 1 

COMMOM /KONST/ Kil) 2 

COMMON /TEMP/ T(U 3 

LUT5=11 ^ 

REWIND LUT5 5 

NN0DES=K(4) 6 

12 READ CLUT5) (T ( I) , 1^= 1 • NNODES ) 7 

RETURN 8 

END 9 


SIBFTC GRAFIC 


SUBROUTINE PICTR 
COMMON/TITLE/ HEAD(20) 

COMMON /TEMP/ T(l) 

COMMON /ARRAY/ IA(1) 

COMMON /KONST/ NSEG, NPAY ,NSURF, NNOOES.MAXNLf NPRI NT»NDAMP,NNL » 

1 NAREAL»NAREANfNARCLN,NTHETA,NITESTf NKHICf NEMISS,NAFA,NSNT, 

2 NEXTRAO), 

3 SIGMAf OPTH, DELS. EMI SS»CONl,CON2. DELDMPtX, TABS, HI KCON,PSWTCH, RAO fH 
DIMENSION NODNO( 15).T0UT(15). LINE(30), NBCD(50), IFMT(13) 
DIMENSION NRENDNUS), NRENDTdS), IFMTK96), Ad) 

EQUIVALENCE (IA,A) 


DATA ILINE.IBLANK /6H******, 
DATA ITARG /6H$SSS$$/ 


6H 


DATA I FMTl ( 1 ) /576H( IHl , 5X. 20 A6/8X. 10 ( IH* ) . 


(8H^ 


) ,2H**, 


C8H ),2H *,10Xf 1H*,15X,IH*/8X 

(8H ).2H *,10X.7H* SPACE. 9X. 

18 ,2H *.10X.10H* NODE NO..I5.2H */ 
)%2H *.10X. 1H*.15X,1H*/8X,10H* TEMPS= 
F8.2 ,2H *,10X,7H* TEMP=,F8.2,2H */8X.10H* EMISIV= , 
F8.4 .2H *.10X.8H* EM IS I V. F7.4. 2H * /8X. 10 ( IH* ) . (8H*** 


110X,17< 1H*)/8X,1H*.9X, 
2,18H* PAYLOAD NODES 
31H* /BXflOH* NODE NOS. 
48X,1H*.9X, (8H 

5 , 

6 


7****») ,2H**.10X,17<1H*)/14HKANGLE OF EACH/14H SEGMENT! DEG* )/ lOX. 

8 F8.1/14HKEMISIVITY OF /14H EACH SEGMENT /lOX. F8,4) 

9 / 

DATA ISU1.ISU2.ITW1.ITW2 
1 / 6HSJRFAC .6HE . 6HTANKWA .6HLL / 

DO 5 I=l.NNODES 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


DATA NBCD( 1)/ 

300H 

1 

2 

3 

4 

5 

6 

7 

8 

14 

1 9 10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

15 

2 20 21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

16 

3 31 32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

17 

4 42 43 

44 

45 

46 

47 

48 

49 

50 

/ 


18 

DATA IFMT(l)/ 

78HC2HP 

8X.A1 

♦30A4/2HP 

A6.A2 

.Al. 

(16 

• IX, 

AD/2 

19 

IHP A6.A2.AI. 

(F7. 

2,A1) ) 

/ 







20 


21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 
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5 T(I)=T(I)+TABS 34 

MAXL= MAXNL+1 35 

NSPACE=NNODES 36 

NWALL=MSPACE-NPAY-1 37 

NCINIT=0 38 

NS1== 1 39 

NS2= MSEC 40 

IFCNSEG .GT. 15) NS2=15 41 

8 NTHl = NTHETA+NSl 42 

NTH2 = NTHETA+NS2 43 

NP1=NSPACE-NPAY 44 

NP2=NSPACE*1 45 

NET1=MEMISS+NS1 46 

NET2=MEMISS+NS2 47 

NEPl=f^ETl+NSEG 48 

NEP2=MEP1+NPAY-1 49 

NESPAC=NEMISS+NSEG+NPAY^1 50 

IFMTU 5)= NBCD(NPAY) 51 

IFMT1(13)= NBCD(NPAY) 52 

IFMT1(24)= N8CD(NPAY-1) 53 

IFMTK35)* NBCD(NPAY) 54 

IFMT1(44)= NBCD(NPAY) 55 

IFMT1(54)= NBCD(NPAY) 56 

IFMT1(64)= NBCD(NPAY) 57 

IFMTl(73)= NBCO(NPAY) 58 

IFMT1C86)= NBCD(NSEG) 59 

IFMT1C95)= NBCD(NSEG) 60 

WRITE(6,IFMT1)HEA0.( J,J=NP1,NP2).NSPACE»(T(I )t I=NP1,NP2),TCNSPACE) 61 

1 »(A(I),I=NEP1,NEP2),A(NESPAC)» ( A ( I ) * I =NTH1 . NTH2 ) » 62 

2 (A(n,I=NETl,NET2) 63 

WRITE(6,2004) H 64 

WRITE(6,2003) 65 

DO 10 I=1,MAXL 66 

M*MAXL-I^1 67 

IM=M-1 68 

9 DO 11 L=1.15 69 

LINE(L)= IBLANK 70 

LlNEa + 15)= IBLANK 71 

NOONO(L)= 0 72 

11 TOUT(L)= 0.0 73 

NS1P1= NSl-1 74 

NC=NC1NIT 75 

N= 0 76 

ISNSW=1 77 

50 IF(M .EQ. 1) GO TO 12 78 

PM1=M-1 79 

52 PM=NITEST+MM1 80 

ITEST=IA(MM) 81 

IFdTEST .LT. 0) ISNSW=-1 82 

12 DO 20 J=NS1.NS2 83 

N2= 2*(J-NS1P1) 84 

Nl= N2-1 85 

JJ=NNL+J 86 

JN=J-NS1+1 87 

NLAY=IA( JJ)+1 88 

IFINLAY .LT. M) GO TO 19 89 

N= N+1 90 

91 

CHECK FOR TOP SURFACE TARGETING 92 

93 

26 IFCNLAY .GT. Ml GO TO 21 94 

C 95 
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NLAY »M THIS IS A TOP SURFACE, CHECK TO SEE IF TARGETED 96 

97 

24 JJ=NEMISS+J 98 

EMISIV=A(JJ) 99 

IFIEMISIV .LT. 0.1) GO TO 21 100 

LINE(M1)=ITARG 101 

LINE(M2)=ITARG 102 

GO TO 27 103 

104 

NUJJ) M THIS IS AN INTERIOR NODE OR TANKWALL SURFACE. 105 

106 

21 LINE(N1)= ILINE 107 

LINE(N2)- ILINE 108 

27 ISIGN=1 109 

IF(I .LT. MAXL) GO TO 23 110 

NRENON(N)sIBLANK 111 

NRENDTIN)=IBLANK 112 

LENO*IBLANK 113 

LLEND=ILINE 114 

NAM1=ITW1 115 

NAM2=ITW2 116 

NAM3=ISU1 117 

NAM4=ISU2 118 

GO TO 18 119 

23 ISIGN-ISNSW 120 

LEND*ILINE 121 

LLEND=ILINE 122 

NAM1*NBC0( IM) 123 

NAM2=IBLANK 124 

NAM3=I BLANK 125 

NAM4=I BLANK 126 

NRENDN(N)=ILINE 127 

NRENDT(N)»ILINE 128 

18 NN= NC+M 129 

NOONO(JN)- NN«ISIGN 130 

TOUT(JN)=T(NN) 131 

19 NC= NC+NLAY 132 

20 CONTINUE 133 

IF(N .GT. 0) GO TO 22 134 

WRITE(6,2001) 135 

GO TO 10 136 

22 IFMT( 6)= NBCD(N) 137 

IFMT(11)= NBCD(N) 138 

WRITE! 6, IFMT)LLEN0, ( L I NE ( L ) , L=1 , 30 ) , NAMl , NAM2 , LEND, 139 

1 (NOONOIL) .NRENDN(L) ,L=1.N), NAM3.NAM4.LEND, 140 

2 (TOUT(L) ,NRENDT(L),L=1,N) 141 

10 CONTINUE 142 

DO 31 1=1,3 143 

31 LlNEtl )=IL1NE 144 

WRITE! 6,2002) (L INE! I ) , I =1 , 3 ) , NWALL, TINWALL), ! L INE (I ) , 1 = 1 ,3) 145 

IF!NS2 .GE. NSEG) GO TO 40 146 

NS1= NS2+1 147 

NS2= NSEG 148 

NCINIT=NC 149 

GO TO 8 150 

40 DO 41 I=1,NN0DES 151 

41 T! I )=T! I )-TABS 152 

RETURN 153 

2000 FORMAUlHl, 5X, 20A6/4 IX, 17 ! IH* ) , lOX, 17 ! IH* ) / 154 

1 41X,1H*,15X,1H*,10X,1H*,15X,1H* / 155 

2 41X, 9H* PAYLOAD, 7X, IH*, 10X,7H* SPACE, 9X.1H* / 156 

3 41X.10H* NODE N0.,I5,2H *,10X,10H* NODE NO.,15,2H •/ 157 
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4 41X.lH*«15X*lH*.10X.lH««15XilH*/ 158 

5 41X.7H* TEMP=,F8.2,2H *, 10X.7H* TEMP=,F8.2.2H */ 159 

6 41XtlH>»15X.lH*tlOXtlH««15X.lH*/ 160 

7 41Xtl7(lH*»,10X,17<lH*) / 161 

8 14HKA'4GLE OF EACH / 14H SEGMENT ( OEG. ) / 10X,15F8.1) 162 

2001 F0RMATJ2HP /2HP /2HP ) 163 

2002 F0RMAT(2HP /2HP /2HP ,41Xt3A6/2HP ,41Xt IH* f 16X, IH*/ 164 

1 2HP .41X.1H**9H TANKWALL.7X.1H* / 165 

2 2HP ,41X,1H*,10H NODE N0.*tI5f2H */ 166 

3 2HP ,41X,lH*t7H TEMP.= F8.2t2H •/ 167 

4 2HP «41XtlH*tl6XtlH*/2HP »41Xt3A6 I 168 

2003 FORMAT (7H LAYER /6H NO. » 169 

2004 FORMAT I25HKTANK TO PAYLOAD SPACING= F7.4, 4H FT. /I 170 

END 171 


$I8FTC PUNTEM 

SUBROUTINE PCHTINI 1 

COMMON /TEMP/ T(l) 2 

COMMON /FIXCON/ FC(1) 3 

COMMON/KONST/ KID ^ 

C SET, NPSWCH TO NON-ZERO TO INDICATE CARDS BEING PUNCHED 5 

K(31)=N 6 

NN0DES=K(4) 7 

FC( 19)=10000. 8 

MRITE(6>101i 9 

PUNCH 100, NNOOES, (T( II, I=1,NN0DES» 10 

100 F0RMAT(18,9F8.1/(10F8.1)I 11 

101 FORMAT! 36HL TEMPERATURE CARDS BEING PUNCHED) 12 

RETURN 13 

END 1^ 
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APPENDIX H 


^^3 


SYMBOLS USED IN ANALYSIS 

..2 2 

area, ft ; m 

2 2 

incremental area, ft ; m 

2 2 

incremental area in x direction in finite difference equation, ft ; m 

2 2 

incremental area in y direction in finite difference equation, ft ; m 

2 2 

incremental area in z direction in finite difference equation, ft ; m 
major semi-axis of tank plus shield, ft; m 

product of 1 surface area times geometric radiation view factor between 
i^^ and surface, ft^; m^ 

minor semi-axis of tank plus shield, ft; m 

specific heat, Btu/(lb)(‘^R); J/(kg)(K) 

geometric radiation view factor - portion of total energy, either emitted or 
reflected diffusely from surface i, that is intercepted by surface j 

conductor value for heat transfer, Btu/(lb)(°F); J/(sec)(K) 

tank-to-payload spacing, ft; m 

thermal conductivity, Btu/(hr)(ft)(°R); J/(sec)(m)(K) 
payload ring width, ft; m 
slope of normal Nj 

normal to center of an incremental surface 
number of surfaces in radiation enclosure 
number of incremental areas in j^^ payload ring 
heat flow, Btu/hr; J/sec 

incremental heat flows, either across a plane of any segment or into (or out of) 
any surface segment, Btu/hr; J/sec 

internal heat generation rate at each node, Btu/hr; J/sec 

distance between any two incremental surfaces, ft; m 

payload radius, ft; m 

distance from tank center to incremental surface, ft; m 
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Rg distance from payload center to incremental surface, ft; m 

S source term, qot/k, (ft^)(°R)/hr; (m^)(K)/sec 

AS thickness of MLl, ft; m 

T temperature, °R; K 

t time, sec 

3 3 

V volume, ft ; m 

X, y, z coordinate system, point coordinates 
major semi-axis of tank, ft; m 
Ax incremental distance in x direction 

y^. minor semi-axis of tank, ft; m 

Ay incremental distance in y direction, ft; m 

Az incremental distance in z direction, ft; m 

a thermal diffusivity, ft /hr; m /sec 

/3 angle measured from minor axis to segment boundary on tank 

y location angle on payload 

e emissivity 

9 angle on tank to specify segments 

p density, Ib/ft^; kg/m^ 

(p angle between the normal to a surface and the vector R from that surface 

angle from top of tank to Rj 
Subscripts: 

i index or counter 

j index or counter 

k index or counter 

lat lateral 

m, n reference or index 

norm normal 

p refers to payload surface increment 

t refers to tank surface increment 
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Notation: 


— vector 

^ unit vector 

2 02 -2 -2 
V Laplacian operator, + — — + 

9x^ dy^ dz^ 

unit vectors in x, y, z directions 
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Figure 6. - Geometric configuration for determining radiation view factors between tank and payload. 
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